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Abstract. We define the correlation of holes on the triangular lattice under periodic 
boundary conditions and study its asymptotics as the distances between the holes grow 
to infinity. We prove that the joint correlation of an arbitrary collection of lattice-triangular 
holes of even sides satisfies, for large separations between the holes, a Coulomb law and 
a superposition principle that perfectly parallel the laws of two dimensional electrostatics, 
with physical charges corresponding to holes, and their magnitude to the difference between 
the number of right-pointing and left-pointing unit triangles in each hole. 

We detail this parallel by indicating that, as a consequence of our result, the relative 
probabilities of finding a fixed collection of holes at given mutual distances (when sam- 
pling uniformly at random over all unit rhombus tilings of the complement of the holes) 
approaches, for large separations between the holes, the relative probabilities of finding 
the corresponding two dimensional physical system of charges at given mutual distances. 
Physical temperature corresponds to a parameter refining the background triangular lattice. 

We also give an equivalent phrasing of our result in terms of covering surfaces of given 
holonomy. From this perspective, two dimensional electrostatic potential energy arises by 
averaging over all possible discrete geometries of the covering surfaces. 



Introduction 

In [3] we considered the joint correlation uj of symmetrically distributed holes on the 
triangular lattice, a natural extension of the monomer-monomer correlation introduced by 
Fisher and Stephenson in [10] . Under the assumption that one of the holes is a unit lattice 
triangle u on the symmetry axis I and all remaining holes are lattice triangles of side 2, 
oriented so that they point away from the lattice line perpendicular to (. that supports 
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w, we proved in [3] that, asymptotically as the distances between holes are large, u 
satisfies a multiplicative superposition principle that perfectly parallels the superposition 
principle of two dimensional electrostatics (with holes corresponding to electrical charges, 
and charge magnitude given by the difference between the number of up-pointing and 
down-pointing unit triangles in a hole) . Our proof was based on explicit product formulas 
we obtained in [4] for the number of lozenge tilings of two families of lattice regions. 

The correlation we used in [3] was defined by including the holes inside a lattice 
hexagon approaching infinite size so that the holes remain around its center. The presence 
of the boundary of the hexagon distorts the local dimer statistics. More precisely, there 
exists an explicit real valued function / defined on the hexagon (/ is the unique maximum 
of a certain local entropy integral; see [8] [7]) so that in the scaling limit the local statistics 
of dimers at each point inside the hexagon is governed by iis,t, where is the tilt of 
/ at that point and iis,t is the unique invariant Gibbs measure of slope (s, t) (this was 
conjectured by Cohn, Kenyon and Propp in [7] and proved and generalized by Sheffield in 
[21] and Kenyon, Okounkov and Sheffield in [17]). It follows from the latter two papers 
that /io,o is the unique invariant Gibbs measure of maximal entropy. Since the function 
/ has tilt (0, 0) only at the center of the hexagon, the local dimer statistics is distorted 
away from maximal entropy everywhere except at the center^: 




(for more details, see [8]; this does not happen in the case of a square on the square 
lattice considered in [10]). Therefore, an important goal is to define the correlation of 
holes in a different way, via regions that do not distort the dimer statistics, and determine 
whether it still reduces to the superposition principle of electrostatics in the limit as the 
holes grow far apart. Other highly desirable features are to allow general, not necessarily 
symmetric distributions of the holes, as well as holes of arbitrary size. 

It is these goals that we accomplish in this paper. We give a new definition for the 
correlation of holes by including them in a sequence of lattice tori of size approaching in- 
finity (see Section 1). We prove that this new correlation also reduces to the superposition 
principle of electrostatics as the distances between the holes grow to infinity. 

We note that by [21] and [17] the torus measures on doubly periodic bipartite planar 
graphs converge to the maximal entropy invariant Gibbs measure. Thus defining hole 
correlations via limits of tori is consistent with our goal of having maximal entropy dimer 
statistics. 

The main result of this paper. Theorem 1.1, addresses the case of an arbitrary collection 
of lattice-triangular holes of even sides. It is a counterpart of our result [3, Theorem 2.1], 
but more general than that in three important ways: (1) no symmetry of the distribution 
of the holes is assumed here; (2) the triangular holes are allowed to have arbitrary even 
sides, and (3) the multiplicative constant has an explicit conceptual interpretation — it is 
obtained by multiplying a simple explicit constant, depending only on the charges of the 
holes, by the product of the correlations of each hole regarded on its own. 

As a consequence of Theorem 1.1, the parallel to electrostatics described in [3], which 
was restricted there by assumptions of symmetry and charge, becomes now established 
in general — including the presence of a "lattice refinement" parameter which corresponds 
to physical temperature. We detail this in Section 15, where we also give an equiva- 
lent phrasing of our result in terms of covering surfaces of given holonomy. From this 
perspective, two dimensional electrostatics arises by averaging over all possible discrete 
geometries of the covering surfaces. 

Regarded from the point of view of surface models, our result raises naturally the 
question of describing the limit of the average height function — which can be suitably 
defined, even though as we see in Section 15 it is multivalued in general — as the lattice 
spacing approaches zero. It turns out that the right object to consider is, by contrast to 
other dimer surface limit shape results in the literature (see [6] [8] [7]), the un-normalized 
average height function (the normalized one approaches zero in the fine mesh limit), and, 
as we show in [5], the limit surface is a sum of helicoids. This in turn is equivalent to the 
statement that a certain naturally defined discrete field (the 90° rotation of the tilt of 
the average height function) converges in the scaling limit to the electric field; sec Section 
15 for more details. The general surface models arising from dimers on doubly periodic 
bipartite planar graphs are studied in [17], where they are classified according to the 
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behavior of the variance of the height function. 

The defining difference between our results and the above quoted ones is the presence 
of holes. Indeed, the boundaries of the holes form the only boundary present in our case. 
If there are no holes, the effects we are studying in this paper and in [5] are not present. 
Consistent with this is that for instance [17] focuses on the variance of height function 
differences, since the behavior of the average height function is not so interesting in the 
absence of boundary. So we could say that when specialized to the hexagonal lattice with 
unit weights the results of [17] give the variance of the height function differences (the 
field mentioned in the previous paragraph being zero), while the sequel [5] of this paper 
gives the value of the field created by the presence of holes. 

It appears that there are very few rigorously proved results on the asymptotics of the 
correlation of non-zero charge holes in the literature (it is this case that reveals electro- 
statics as the governing law; the particular instance of holes of charge zero equivalent to 
dimer-dimer correlations has been studied more; see e.g. [10]). In fact, the only such re- 
sult the author is aware of is the monomer-monomer correlation on the square lattice along 
a lattice diagonal direction obtained by Hartwig [11] (this was conjectured in [10]). For 
a discussion of the differences between our set-up and other discrete models for Coulomb 
gas in the physics literature, and for some earlier pointers to electrostatics suggested by 
the study of the spin-spin correlation in the Ising model, see [3, pp. 2-3 and 91]. 

In order to relate he results of this paper relate to those of [2] and [3] , we note that 
despite their very different definitions, the correlation Cj considered in this paper and 
the correlation u) of [3] turn out to satisfy, up to a multiplicative constant, the same 
asymptotic superposition principle. We also note that the special case of Theorem 1.1 
when one has just two oppositely oriented lattice triangular holes of side 2 was treated 
in [2] using a third, quite different definition of correlation, and was shown there to obey 
exactly the same asymptotics as the one that follows by Theorem 1.1 for a; in this special 
case. 

We conclude by mentioning that, as pointed out by the anonymous referee, in a differ- 
ent surface model, the discrete Gaussian free field model (see e.g. [22]), there is a natural 
analog of holes. Moreover, the question of determining the scaling limit of the average 
height function in the presence of these "holes" is much simpler than what it turns out 
to be in our model, and also leads to two dimensional electrostatics. In view of previous 
connections between the discrete Gaussian free field model and the dimer model (see e.g. 
Boutillier [1] or Kenyon [16]), this may offer an intuitive reason why the results we obtain 
in this paper and in [5] should hold. 

1. Definition of Co and statement of main result 

Draw the triangular lattice so that some of the lattice lines are vertical. For terminolog- 
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ical brevity, we will often refer to unit lattice triangles as monomers; a left-pointing (resp., 
right-pointing) unit lattice triangle is called a left-monomer (resp., right-monomer) . 

For any finite union Q of unit holes on the triangular lattice, define the charge ch.{Q) 
of Q to be the number of right- monomers in Q minus the number of left- monomers in Q. 

Let Ui,...,Un be arbitrary unit holes on the triangular lattice. Following Kenyon 
[14], provided Yl7=i ch(?7i) = 0, we define the joint correlation of t/i, . . . , t/^ by 

fTT TT \ V M(gAr,Ar \ ^1 U ■ ■ ■ U ^^) 
c.,(C/i, ...,[/.):= Jim ^ , (1.1) 

where H^ n is a large lattice rhombus of side N whose opposite sides are identified so as 
to create a torus, and M{R) denotes the number of lozenge tilings of the lattice region 
R (a lozenge, or unit rhombus, is the union of any two unit lattice triangles that share 
an edge)^. 

An important advantage of this definition is that, provided Ui U ■ ■ ■ U Un is a union 
of non-overlapping lattice triangles of side 2, [14] and [15] provide an expression for 
u>i{Qi, . . . , Qn) as a determinant. 

To state this explicitly, we introduce the following system of coordinates on the mo- 
nomers of the triangular lattice. Choose the origin O to be at the center of a vertical 
unit lattice segment, and pick the x- and y-coordinate axes to be straight lines through 
O of polar arguments —tt/3 and tt/3, respectively. Coordinatize left- monomers by the 
(integer) coordinates of the midpoints of their vertical sides; coordinatize each right- 
monomer likewise, by the coordinates of the midpoint of its vertical side (thus any pair of 
integers specifies a unique left-monomer and a unique right-monomer, sharing a vertical 
side) . 

Let (liJ'i), . . . , (/fe, I'lJ be the coordinates of the left-monomers contained in the union 
t/i U • • • Ut/^i of our holes, and let the coordinates of the right- monomers contained in this 
union be {ri,r[), . . . , (r/,, r^) (the latter are the same in number as the former, since we 
are assuming the total charge of the unit holes to be zero). Then, provided UiU ■ ■ - UUn 
is a union of non-overlapping lattice triangles of side 2, it follows by [14] and [15] that 

ujiiUu ...,Un) = \det (Pin - lj,r[ - lj))i<i,j<kl (1-2) 

where 



^TT^ Jo Jo l + e-^^ + e-*-^ ^ ^ 



■^The total charge of the the holes needs to be zero in order for the region in the numerator of (1.1) 
to contain the same number of each type of unit triangles — a necessary condition for the existence of 
lozenge tilings. For any situation when the total charge is not zero, (1.1) would assign value to the 
joint correlation of the holes, irrespective of their relative positions. 
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Thus, if Qi, . . ■ , Qn are holes that are unions of non-overlapping triangles of side 2, 
i = 1, . . . , n, and Y17=i ch(<5i) = 0, defining their correlation by coi seems both natural 
and convenient, as the above determinant formula can be used as a starting point for 
determining its asymptotics. 

What if Yl^=i ch(Qi) 7^ 0? We need then an appropriate extension of the correlation 
uji that applies also in this situation, and is manageable enough to work with so that 
its asymptotics can still be worked out. We found that the following inductively defined 
correlation a) does the job. This was a crucial part of our proof. 

Definition. Let Ti, . . . , be lattice triangular holes of side 2. 

If q — ^iLi ch(Ti) > 0, define their joint correlation co{Ti, . . . ,T„) inductively by: 
(i) If q — 0, define Cj by (1.1). 
(u) If q = 2s, s > 1, define 

ob{Ti,..., T„) = lim (3r)2^cD(Ti, . . . , T„ W{3r, 0)), 

1 — >oo 

where W{3r, 0) is the left- {or west-) pointing lattice triangle of side two whose western- 
most left-monomer has coordinates (3r, 0) . 

If q < 0, consider the mirror images Ti, . . . ,Tn of our holes across a vertical lattice 
line, and define a>(Ti, . . . , T„) := a>(T'i, . . . , T^). 

(The fact that u) is well-defined by the above follows from Propositions 3.2 and 4.1.) 

Note that lj is defined via a sequence of toroidal regions of sizes growing to infin- 
ity. Since, in stark contrast with the case of the hexagonal regions used to define the 
correlation co of [3], the periodic boundary conditions of such regions do not distort in 
the limit the local dimer statistics, the definition u for the correlation seems much more 
natural from a physical point of view than the u defined in [3]. The question is: Does the 
asymptotics of uj also lead to the superposition principle of electrostatics? The object of 
this paper is to prove that the answer to this question is affirmative. 

To state our results explicitly, pick a reference monomer /j, in each hole Q under 
consideration, and define Q{a, b) to be the translated copy of Q in which the image of jj, 
has coordinates (a, 6). 

The main result of this paper is the following. 

Theorem 1.1. Let Ai, . . . , be lattice-triangular holes of even side-lengths {see Figure 
1.1(a) for an illustration). Then for any family of n distinct pairs of integer multiples 
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of 3, (xi, Hi), i = 1, . . . ,n, we have 



/ X -HEr=i I ch(Aoi-i Er=i ch(A,)i} n 

X n d(A,(i?a;,, , A.' (^^i , ^^.))| ch(A,)ch(A,) 

l<i<j<n 

+ 0{R^ Ei<i<,<„ ch(Ai) ch(A,))-l>j^ 

as R oo, where d denotes the Euclidean distance. 

To prove this, it turns out to be convenient to generalize first the above statement by 
introducing a mnnber of new parameters. These will be essential for the evaluation of a 
certain determinant which represents a key step of our calculations. 

To state this generalization, let Qi be an arbitrary union of non-overlapping lattice 
triangles of side two, for z = 1, . . . , n. Since Qi is not necessarily connected, we call it a 
multihole. 

Consider the case of multiholes satisfying the following three conditions: (1) each Qi is 
pure, i.e., is a union of like-oriented lattice triangular holes of side 2; (2) each Qi is linear, 
i.e., its constituent lattice triangular holes of side 2 have their centers along a straight 
(not necessarily lattice) line; and (3) the directions of these pure and linear multiholes 
are all parallel among themselves (an instance of such multiholes is pictured in Figure 
1.1(b)). 

We will say that a rational number r is divisible by 3 if r = a/6, where a and b 
are relatively prime integers and 3|a. The generalization of Theorem 1.1 we need is the 
following. 

Theorem 1.2. If Li, . . . , are pure, linear and parallel multiholes, and if the slope 
q & Q of their common direction satisfies 3\1 — q, then for any family of n distinct pairs 
of integer multiples of 3, {xi,yi), i = 1, . . . ,n, we have 

HEr=i|ch(LO|-|Er=iCh(LO|} n 

'] 

27T 



Cj{Li{Rxi,Ryi),...,Ln{Rxn,Ryn)) ^ [-^] W'^i^i) 



i=l 



X J] d{Li{Rxi, Ryi), L,(i?a;,-, i?y,))^ ^'^(^^) 

l<i<j<n 



as R ^ oo, where d denotes the Euclidean distance. 
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(a). (b). 

Figure 1.1. (a) Triangular holes of even side, (b) An example of pure, linear, parallel 
multiholes (ellipses indicate that multiholes keep fixed size and shape as they 

move away from one another). 



We note that the assumption 3|1 — g on the slope q of the common direction of the 
linear multiholes is not really essential, but makes the calculations in Sections 5-7 simpler 
(outside this assumption, Newton's divided differences need to be adjusted by introducing 
some multiplicative constants). In addition, the set of rationals q with 3|1 — g is easily 
seen to be dense in Q. Similarly, the divisibility by 3 of the integers Xi and j/j is not 
really essential, but simplifies the application of the method of factor exhaustion to the 
evaluation of the determinant M" given by (5.11)-(5.15) (sec Sections 9, 11 and 12). It 
clearly does not impend on the generality of the directions of the rays through the origin 
that contain (xj,|/j). 

Geometrically, the above result states that if the multiholes Li, . . . , L„ are translated 
away from each other at the same rate along fixed rays through the origin, with their 

shapes and sizes kept intact, then the asymptotics of their joint correlation u is given 
by an explicit multiplicative superposition principle. This corresponds to a refinement 
of the superposition principle of electrostatics: the -R-dependent part on the right hand 
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Figure 1.2. An instance of pure, linear multiholes of arbitrary directions. 

side of (1.4) parallels perfectly its physical counterpart (our statistic ch corresponding to 
electrical charge), but our multiplicative constant is given by the product of the correla- 
tions of the constituents of each multihole (this would correspond to the physical charges 
having some "structure" ) times an absolute constant raised to a power that depends just 

on the ch(Li)'s. 

Equivalently, we could hold the location of the multiholes fixed, and refine the triangu- 
lar lattice instead (more precisely, for each i, fix the geometrical position of some lattice 
point Pi in Lj and shrink Lj around Pi in step with the refinement of the lattice). From 
this perspective. Theorem 1.2 states that the scaling limit of the joint correlation of pure, 
linear and parallel multiholes is given by the described refined superposition principle. 
As the lattice spacing approaches zero, the multiholes, preserving their charge, become 
more and more point-like. Nevertheless, the specific way in which each multihole Lj 
approaches this point-likeness — the structure of Lj — makes itself felt as a multiplicative 
contribution of cD(Lj). 

Theorem 1.1 follows easily from Theorem 1.2. 

Proof of Theorem 1.1. By Lemma 3.1(b), due to forced unit rhombic tiles, a triangular 
hole of side 2s is equivalent to a horizontal string of s contiguous triangular holes of side 
2 (an instance of this is illustrated in Figure 1.1(b): the union of the leftmost three holes 
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in Qi forces a total of 12 unit rhombic tiles; their removal creates a triangular hole of side 
6). Since in our system of coordinates horizontal lines have slope q = 1, the statement 
follows by Theorem 1.2. □ 

Although the statement of Theorem 1.2 is general enough so that we can prove the 
determinant evaluation it boils down to by factor exhaustion (and thus deduce Theorem 
1.1), it is worth mentioning that we can extend it further to the case when the pure, 
linear multiholes have arbitrary, not necessarily parallel directions. 

Theorem 1.3. The statement of Theorem 1.2 holds more generally when the slopes q^ 
of Li and q'j of Rj are arbitrary rational numbers with 3\1 — q^, 3\1 — q'j, i = 1, . . . ,m, 
j = l,...,n. 

Such a more general instance is pictured in Figure 1.2. 

Remark 1.4. The approach employed in this paper does not allow the presence of holes 
shaped as lattice triangles of odd size. This is due to the fact that a certain advantageous 
employment (noticed by Kenyon [15]) of a construction of Percus [20] for obtaining a 
Pfaffian orientation of a subgraph from a Pfaffian orientation of the full graph, which 
allows expressing the ratio between the number of perfect matchings of a subgraph and 
the number of perfect matchings of the full graph as a determinant, is not applicable in 
this situation. 

2. Deducing Theorem 1.2 from Theorem 2.1 and Proposition 2.2 

Denote by E{x, y) the right- (or east-) pointing lattice triangular hole of side 2 whose 
central monomer (a left-monomer) has coordinates {x,y). For any g e Q and any list 
a = [ai, . . . , an] of integers for which qai e Z, and the E{ai, qaiYs are mutually disjoint, 
z = 1, . . . , n, define the multihole E^ by 

El := E{ai, qai) U • • • U E{an, qan). 

Similarly, let W{x,y) be the left- (or west-) pointing lattice triangular hole of side 2 
whose central monomer (a right-monomer) has coordinates (x^y). Given g G Q and a 
list of integers b = [6i, . . . , 6„] so that qbi G Z, and the W{bi, qbiYs are mutually disjoint, 
z = 1, . . . , n, define the multihole by 

:= W{buqbi) U • • • U W{bn, qbn). 

By definition, any pure and linear multihole is of the form E^ or W^. The parameter 
q gives the slope of the line containing the centers of the constituent triangular holes of 
side 2. Thus, the members of any family of pure, linear and parallel multiholes have one 
of the above two forms, and they all share the same q. 
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Theorem 2.1. Let m,n >0 be integers, and q a rational number for which 3|1 — g. For 
i = 1, . . . , m, let SLi = [an < ■ ■ ■ < ais^] be a list of Si integers so that qaij G Z and the 
E{aij, qaij) 's are mutually disjoint, j — 1, . . . ,Si. For i — 1, ... ^n, let hi = [bn < ■ ■ ■ < 
bit^] be a list of ti integers so that qbij e Z and the W{bij, qbij) 's are mutually disjoint, 

Then for any family of m + n distinct pairs of integer multiples of 3, {xi,yi), i = 
1, . . . ,m, (zjjWj), j = 1, . . . ,n, we have, as R ^ oo, that 

u{El^{Rxu Ryi), EljRxm,Rym), W^Rzi, Rwi), W^jRZn, RWn)) = 
/-\ E^i *i+Er=i EIli *i-Er=i til 

m n 

^ n n (^^i ~ ^ikf n n ^^^^^ ~ ^^'^)^ 

i=l l<j<k<Si i=l l<j<k<ti 

^ 

YlZi U^^iii^i - ZjY + [xi - Zj){yi - Wj) + {yi - w^YY^^i 

X n [{Zi-Zjf + {Zi-Zj){Wi-Wj) + {Wi-Wjf]'''i 

^2{El<i<j<m "l"El<i<j<T!, *i*J~"Ei = l Ej = l 

Proposition 2.2. Let q e Q with 3|1 — g, and let a = [ai, . . . , as] be integers so that 
qai e Z and the E{ai, qai) 's are mutually disjoint, i = 1, . . . , s. Then 

'^(^I)=f||) + q + n K (2-2) 

\ / l<i<i<s 

Theorem 1.2 follows easily from the above two results. 

Proof of Theorem 1.2. Re-denote the integer n in the statement of the theorem by 
N. Assume that m of the given L^'s point east; by assumption, they are of the 
form £"1^, . . . , -E"!^, where g G Q, 3|1 — g, and = [aa, • • • , ftisj are lists of integers. 
The remaining n := N — m multiholes are then of the form W^^, . . . , W^^, where = 
[bii, . . . , biti] are lists of integers. 

Then the left hand side of (1.4) is clearly the same as the left hand side of (2.1). By 
Theorem 2.1, to complete the proof it suffices to show that the right hand side of (2.1) is 
equal to the right hand side of (1.4). 

11 



Clearly, by our definition of cD in Section 1, Proposition 2.2 implies that for any q E 
3|1 — g, and any list b = [bi, . . .,bt\ of integers for which is defined, one also has 



^(n)=(^] (i+q+q'P n (^^-^.)'- (2-3) 

\ / ^<i<j<t 

By (2.2) and (2.3), the expression on the right hand side of (2.1) can be rewritten as 



2ti 



ni<i<,-<m[(^^^ - ^^j)^ + i^^i - R^j){Ryi - Rvj) + {Ryi - Ry. 



1 ]Tj=i[{R^i - Rzj? + {R^i - Rzj){Ryi - Rwj) + {Ryi - Rwj)^Y^'^ 
X Yl [(R^i - R^if + iR^i - Rzj)iRwi - Rwj) + (Rwi - RwjfY'*^ . 

l<i<j<n 

However, since in our 60°-angle system of coordinates the Euclidean distance is given by 
the formula 

d((a, b), (c, d)) = ^/{a - cf + (a - c)(6 - d) + (6 - df, 

the quantities in the square brackets are recognized to be the squares of Euclidean dis- 
tances between the reference monomers in the pairs of multiholes 

{ElXRx,,Ry,),El.{Rxj,Ryj)\ 
[yVl {Rz„ Rw,), {W^^ {Rzj, Rwj)), 

and 

{El {Rxi, Ryi), < . {Rzj, Rwj)), 

respectively. 

Furthermore, it is clear by definition that ch(£^|.) = 2si and ch(W^.) = —2tj, for 
i — l,...,m, j = l,...,n. Therefore the expression above is just the right hand 
side of (1.4) when the collection of multiholes {Li,...,Liv} is written out explicitly 
as {E^, . . . , -E"!^^, W^bi' • • • ' }• This completes the proof of Theorem 1.2. □ 

The above arguments reduce the proof of Theorem 1.2 to proving Theorem 2.1 and 
Proposition 2.2. 
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Outline of proof of Theorem 2.1. The arguments in the proof can be grouped in three 
parts. 

The first part, covered by Sections 3 and 4, proceeds by providing an expression for 
Cj (see Proposition 3.2) as the determinant of a certain matrix M whose entries are of 
two types: on the one hand, values of the function P defined by (1.3), and on the other, 
values of certain two-variable polynomials Us{a, b), s > 0, that are the coefficients of the 
asymptotic series of P{—R — 1 + a, — 1 + 6), as R ^ oo. Proposition 4.1 provides an 
explicit expression for the Us{a,b)^s in terms of powers of the finite difference operator. 
We obtain this expression by Laplace's method for finding the asymptotic series of a 
contour integral. We use in our derivation Gauss' summation of a 2F1 hypergeometric 
function and Newton's expression for a polynomial in terms of powers of the difference 
operator. 

To prove the statement of Theorem 2.1 one needs then to determine the large R 
asymptotics of det(M). However, the matrix formed by the leading i?-parts of the entries 
of M turns out to be singular. The second part of our proof, contained in Sections 5-7, 
addresses this complication by finding explicit elementary row and column operations on 
the matrix M so that the resulting matrix M' — clearly of the same determinant as M — 
has the property that the matrix formed by the asymptotics of its entries is non-singular. 
The row and column operations that achieve this turn out to be governed by Newton's 
divided difference operator. It is in this part that we crucially use the pure, linear and 
parallel structure of our multiholes. 

Finally, in the third part (covered by Sections 8-12) we evaluate the determinant of 
the matrix M" formed by the asymptotics of the entries of M' by the method of factor 
exhaustion. This well-known method of evaluating determinants with polynomial entries 
applies when the value of the determinant is guessed explicitly, and is a polynomial with 
all roots and multiplicities known. The method of factor exhaustion proceeds by show- 
ing that the determinant admits all these guessed roots with the guessed multiplicities, 
by finding a suitable number of independent row or column combinations that vanish 
when the variable equals the root in question. This, together with a comparison of the 
degrees, implies that the value of the determinant is equal to the guessed polynomial up 
to a multiplicative constant. An additional argument is then needed to establish this 
multiplicative constant, and thus the value of the determinant. 

The guessed value of det(M") follows from the expression on the right hand side of 
(2.1). The entries of M" turn out to be rational functions in the variables q and Xj, i/i, 
Zj, Wj, i = 1, . . . ,m, j = 1, . . . ,n. In Section 9 we prove that, regarded as a polynomial 
in g, det(M") is divisible by the factors involving q in its guessed expression. Comparing 
degrees in q we deduce that it suffices to prove the determinant evaluation for the matrix 
Mq obtained from M" by specializing q = 0. To this end, we apply the method of factor 
exhaustion to dMo, where d is the least common multiple of the denominators of the 
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entries of Mq. 

We provide vanishing row and column combinations that prove the divisibihty of 
detdMo by the factors of its guessed expression in Sections 11 and 12. The fact that 
the specified row and column combinations are vanishing is deduced from known exact 
evaluations of 2-^1 and 3F2 hypergeometric functions. This allows us to deduce in Section 
10 the validity of the determinant evaluation up to a multiplicative constant. Finally, in 
the Appendix we show that the multiplicative constant is the one dictated by the right 
hand side of (2.1). 

As described above, the obtained evaluation of det(M") determines the asymptotics 
of det(M), and thus that of a). This leads to the expression on the right hand side of 
(2.1), and completes the proof of Theorem 2.1. The details are given in Section 13. □ 

Outline of proof of Proposition 2.2. We use the approach described in the first two 
parts of the above outline of the proof of Theorem 2.1, applying it to the single multihole 
El^. The corresponding matrix M" turns out to be an upper-triangular block-matrix 
consisting of 2 x 2 blocks. The determinants of the diagonal blocks are readily calculated, 
and their product is checked to agree with the expression (1.6). The detailed proof is 
presented in Section 13. □ 

The necessary changes in the proof of Theorem 1.2 that afford the more general The- 
orem 1.3 are presented in Section 14. 



Recall the toroidal regions Hjq,N considered at the beginning of Section 1. We say that 
a collection lA = {ui, . . . , Un} of holes on the triangular lattice forces a unit rhombic tile 
t on the lattice if t is present in all tilings of Hn^n \ {-ui, . . . , w^}, for all large enough 
that Hn^n contains U. In proving our determinant formula for a) we make use of the 
following simple observation. This is also needed to deduce Theorem 1.1 from Theorem 



Lemma 3.1. (a). Let U = {ui, . . . ,Un} be a collection of unit holes having total charge 
zero, and let ti, . . . ,tk be unit rhombic tiles that are forced by U. Then 



(b). // the collection {Ti, . . . , T„} of non- overlapping lattice triangular holes of side 2 
forces the unit rhombic tiles ti, . . . ,tk, then 



3. A determinant formula for a) 



1.2. 



Wi(Wi, . . . , Un, tl, . . . , tk) = Wl(Wl, • • • , 



Un)- 



a)(Ti, . . . , Tn, ti,..., tk) — a)(Ti 




Proof. Both statements follow directly from the definitions of uoi and uu in Section 1. □ 
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Consider the function P defined by (1.3). For a,b & Z and s > 0, define Us{a,h) by 
writing the asymptotic series ofP(— 3r — 1 + a,— l + 6),r— >oo, as 



P(-3r - 1 + a, -1 + 5) ~ ^(3r)-^-it/,(a, h). 

s=0 



By the definition of asymptotic series (see e.g., [19, p. 16]), this means 



n-l 



P(-3r - 1 + a, -1 + 6) - ^(3r)-"-^C/^(a, h) = 0{i 



(3.1) 



s=0 



for each fixed value of n. The fact that P(— 3r — 1 + a, — 1 + 6) admits an asymptotic 
series follows from Proposition 4.1, which also gives an explicit formula for the ?7s(a, 6)'s. 

For notational brevity, it will be convenient to define the following families of 2 x 2 
matrices: 

^P(x-l,y-l) P{x-2,y) 
P{x,y-2) P{x-l,y-l) 



A{x,y) :-- 



(3.2) 



and, for s > 0, 



Bs{x,y) : = 



Us{x,y) Us{x-l,y+l) 
Usix+l,y-l) Usix,y) 



(3.3) 



The following result extends the determinant expression (1.2) for coi to a). In particular, 
it provides a determinant formula for the exact correlation of an arbitrary collection of 
non-overlapping triangular holes of side two ((1.2) addressed this only when there were 
the same number of holes of each orientation in the collection) . 

As in the previous section, let E{x,y) and W{x,y) denote the lattice triangular holes 
of side 2 that point east and west, respectively, and have their central monomer of coor- 
dinates {x,y). 



Proposition 3.2. Form > n we have 



), W{ci, di),..., W{cn, dn)) = |det M 
15 



where M is the m x m matrix 

A{ai—ci,bi-di) A{ai—C2,bi—d2) ■■■ A(ai-c„,6i— d„) Bo{ai,bi) Si(ai,&i) ••• Bm-n-i{ai,bi) 
A(a2-ci,b2-di) A{a2-C2,b2-d2) ■■■ A(o2-c„,62-d„) Bo(o2,62) -61(02,62) ••• Bm-n-iiaiM) 



- A{am—Cl,bm—dl) A{am—C2,bm—d2) ■■■ A{am — Cn,bm—dn) Bo{am,bm.) Bi{am,bm) •■• Bm-n-l{am,bm) 

(3.4) 



Proof. We proceed by induction on m — n. Let m — n = 0. By definition, when the 
total charge of the holes is zero, lj is the same as ivi. Therefore, by Lemma 3.1(a), the 
value of tD on the left hand side above remains unchanged if we replace each E{ai, bi) by 
the pair of right-monomers of coordinates (a^ — 1, 6^) and (a^, 6j — 1), and each W{cj, dj) 
by the pair of left- monomers of coordinates (cj,6j + 1) and {cj + l,6j). Since in this 
collection the 4m monomers come in adjacent pairs. Theorem 2.3 of [15] is applicable, 
and, combined with (1.2), provides an expression for their joint correlation Co as 



det (P(r,-/,.,r^-/;.))i 



<i<j<2m 



where (r^, and {k, ZQ, z = 1, . . . , 2m, are the coordinates of the right-, respectively left- 
monomers in the collection. It is immediate to check that when listing the coordinates 
of the right-monomers in the order [(a^ — 1, bi), {tti, 6^ — 1) : i = 1, . . . , m], and those of 
the left-monomers in the order [(c,, bi + 1), (c^ -|- 1, 6^) : i = 1, . . . , m], the above matrix 
becomes precisely the m = n specialization of matrix (3.4). 

Suppose now the statement holds ior m — n — k. By the definition of tD and by the 
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induction hypothesis we have 

Cj{E{ai, 61), . . . , E{an+k+l, &n+fc+l), W{ci,di), . . . , W{Cn, dn)) 

= hm {Srf'+^uiEiai, 61), ... , E{an+k+i, &n+fe+i), W{ci, di), W{cn, dn), W{3r, 0)) 

r— >oo 

= hm (3r)2'=+2 

1 — >oo 

A(ai-ci,6i-di) ••• A{ai-Cn,bi-dn) A(ai-3r,6i) Bo(ai,6i) ••• Bfc_i(ai,6i) 



X 



det 



A{a2-ci,b2-di) 



A{aN—ci,bM—di) 



A{a2-Cn,b2-dn) A{a2-3r,b2) Bo(a2,62) 



A(aN — Cn,bN—dn) j4(oiv— 3r,6jv) -Bo(ajv,6jv) 



-Bte_i(a2,62) 



Bk-i{aN,bN) - 



(3.5) 



where = n + /c + 1. By (3.1)-(3.3), one sees that adding suitable multiples of the last 
2k columns to columns 2n + 1 and 2n + 2, the 2x2 block at the intersection of rows 
2i — 1, 2i and columns 2n + 1, 2n + 2 in the above matrix becomes: 



(3r)-^-iC/fc(ai, h) + 0(r-^-2) {Zr)-^-^Uk{ai - 1, 6^ + 1) + 0(r-^-2) 

(3r)-^-iC/fc(ai + 1, 6i - 1) + 0{r-^-'') {Zr)-^-^Uk{a,, h) + 0(r-^-2) 

(3.6) 

for i — 1, . . . , N. The contributions of the terms in the expansion of the determinant 
in which at least one of the 0{r~'^~^) parts is chosen is clearly 0{r~^^~^). Therefore, 
omitting the 0{r~^~'^) terms in these two columns does not change the limit (3.5). After 
omitting these terms, the factors {3r)~''~^ can be factored out along these two columns, 
cancehng the factor (3r)^'^+^ in the limit (3.5). By (3.3), the 2x2 block resulting this 
way from (3.6) is just Bk{ai, hi). Interchange the pairs of columns {2n + 1, 2?i + 2/c + 1) 
and {2n + 2, 2?i + 2/c + 2) of the matrix resulting this way from the matrix in (3.5) to 
complete the induction step. □ 
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4. An exact formula for Us{a,b) 



The goal of this section is to prove the following result. 

Proposition 4.1. Let P be given by (1.3), and let a,b eZ. Then the coefficients of the 
asymptotic series 



P(-3r - 1 + a, -1 + 6) ~ ^{3r)-'-^Us{a, b), r ^ oo 



s=0 



are given by 



Us{a,b) = -^ [C-'-\l - DC')-' - C''+'^\l - DC)-'] (4.1) 

where D is the difference operator {for a function f defined on Z, Df is defined by 
Df{x) = fix + 1) - fix)), and ( = e^'^*/^. 



Consider the functions 



p(t) = -ln(-l-t) 
Q(t) =t-^(-l-t)-«. 



(4.2) 

(4.3) 



In proving Proposition 4.1 we will use some properties of these functions, which are 
stated in the following three lemmas. 



Lemma 4.2. The power series expansion of Qit) around t = z is 

k 



Qit) = z-\-l-z)-J2 



k>0 



(-1)^ 

klzi'i 



'^j:{'l)ia + b + l)k-l{b)iz^-^ 



{t-z)^. (4.4) 



Proof. By Taylor's formula, it suffices to show that 

k k 



^t-\-l-t)-'' 
dt^ ^ ' 



= z-^-l-z)-" 



(4.5) 



One readily sees that 



s 



t-' = i-iyib)it 
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-b-l 



and 



{-l-t)- = {a)i{-l-t) 



—a— I 



for / > 0. Therefore, by Leibniz's formula for the derivative of a product, we obtain 



t=z 1=0 



Thus, to prove (4.5) it suffices to show that 



E (;) (^)^ ^"'(1 + ^)"^'"'^ = (1 + ^)"' E (;) + ^ + 0/^-^ (4-6) 

Multiplying by (1 + ^)'^, both sides above become polynomials in z~'^ of degree at most 
k. Therefore (4.6) will follow provided we check that the coefficient of z~'^ is the same 
on both sides, for all < n < A;. This is readily seen to amount to showing that 



E 

l=n 



k\ ( I 



1 1 \ l — n 



{a)k-i {h)i = (a + 6 + n)k-n {b)r 



(4.7) 



for all < n < fc. 

The left hand side of (4.7) can be expressed in terms of a hypergeometric function^ by 
writing 

k\ ^ (-l)'(-fc); ^ (-l)'(-fc)^ (-fc + n);-„ 
I U 



I 



n\ {n+l)i-n 
n\{n + l)i-n (n + l)z_ 



l — nj n\{l — n)\ n\{l — n)\ {I — n)\ 



^The hypergeometric function of parameters a\, . . . ,ap and b\, . . . ,bq is defined by 

ai , . . . , ap 



pFq 



hi,..., hq ' ' 



^ (ai)fc ■ ■ ■ (ap)fc __fc 

^^k\{hl)k---{hq)k ' 



where (a)o := 1 and (a)fe := a(a + 1) ■ ■ ■ (a + fc — 1) for fc > 1. 
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ia)k-i (a)fe- 



(a)k-i 



{a)k-n {-a - k -\- n -\- 
{h)i = ib)nib + n)i-n- 
Indeed, the left hand side of (4.7) becomes 

k 



E 



I I \l — n 



{a)k-i {b)i 



{-ir{-k)n{b)n{a)k-n 



{-ir{-k)n{a)k-n{b), 



E 

l=n 



{-k + n)/_„ {n + (6 + n);-^ 
(n + l)i-n {I -n)l {-a - k + n+ l)i-n 



—k + n,b + n 
-a — k + n + 1 



(4.8) 



By Gauss' formula (see e.g. [23, (1.7.6), Appendix (III. 3)]), for a, 6, c e C with Re (c 
a — b)>0 and c 7^ 0, —1, —2, . . . , one has 



a, b 



r(c)r(c-a-6) 
r(c-a)r(c-6)' 



(4.9) 



Let b, k > n he fixed, and choose a to be a large enough negative integer so that 
-a-b-n+l>0 and -a-k + n+ l>0. Then by (4.9) we get 



2-fl 



—k + n,b + n 
-a — k + n + 1 



r(-a +l-{k-n)) T{-a - b - k + 1 + {k - n)) 



r(-a+l) 
1 



(-a+ 1 - A; + n)fc_n 



r(-a-6-A; + l) 
{-a-b- k + l)k-n, 



where at the last equality we used repeatedly that r{x + 1) = xT^x). Therefore, by (4.8) 
equality (4.7) amounts to 



(-l)-(-/c)n,,, , , {-a-b-k + l)k- 
-{b)n {a)k- 



k\ 



{-a + l-k + n)k-n 



n 



{a + b + n)k-n {b)n- 



Dividing by (;) 



and then multiplying by the denominator leads to having 



the same factors on both sides. This proves (4.7) for large enough negative integers a. 
Since both sides of (4.7) are polynomials in a, it follows that (4.7) holds for general a. 
This completes the proof of the Lemma. □ 
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Lemma 4.3. For z ^ —1 the function p{t) given by (4.2) is analytic in a neighborhood 
of z, and the reversion of its power series expansion 



is given by 



V := p{t) - p{z) =pQ{t- z)^- pi{t - zY + 



t-z = {1 + z){e-'' -1). 



(4.10) 



Proof. The first part of the statement is clear. The expression for the reversion of the 
series also follows easily by exponentiating —v = —p{t) +p{z) = ln(— 1 —t) — ln(— 1 — z): 



-1-t 
-1-z 



e"^ ^ t - 2 + 2 + 1 = (1 + 2)e-^ ^t-z = {l + z){e-'' - 1). 



Clearly, for the function p{t) given by (4.2), p'(t) = — 1/(1 + 1), so 

1 -(l + z)-(t-^). 



p'{t) 



□ 



(4.11) 



Lemma 4.4. Replacing t — z in (4.4) and (4.11) by the expression given by (4.10) leads 
to 

Q{t) 



^ = z-\-i-z)-y 



1=0 



x=a+h—l 



(4.12) 



where D is the difference operator. 

Proof. Denoting by [v^]f{v) the coefficient of in the series /(v), we need to show 
that for all A; > one has 



p\t) 



s>0 \.l=Q 



{h)iz 



(1-e-T 



X 



k 

{-(1 + z)-{l + z){e-- - 1)} = -z-\-l - ^)-'^^ E ^^D\x^) 

1=0 



x=a+b—l 



(4.13) 



Since 



(-j-l)V {-!)%. ,,fe 



t\ 



k\ 
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and thus 



i=0 ^•'^ 3=0 

we obtain that the left hand side of (4.13) equals 



s>0 



3=0 



X 



p^(^'j^{a + b + l)s-i {h)iz- 



(4.14) 



Thus, equality (4.13) is equivalent to 



s>0 



E(-i)' f •) + 1)' E f / + ^ + 0.-^ (^)^ ^" 

i=0 11=0 



1=0 



.-I 



x=a+b—l 



which in turn, by extracting the coefficients of 2; and dividing by {—l)^{h)m amounts 
to 



s>0 



3=0 ^-^^ 



n 



(a + 6 + n),_„ = i^"(^'=) ' (4-15) 



for all n > 0. Using the formula for the powers of the difference operator (see e.g. 
[24,p.36]) 



^v(«) = E(-ir^(')/(«+i), 
3=0 ^•'^ 



(4.16) 



the sum in the brackets in (4.15) is readily seen to be just {—lYD^x^\x=i. Therefore 
(4.15) becomes 



E 

s>n 



{-ly fs 



(-1)' 



{a + h + n)s-nD^x''\^=i = 
s\ \nj nl 
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n^k I 



|a;=a+6— 1 • 



We show more generally that for any polynomial f{x) one has 

J2^-^(')ic+ri)s-nD'f{x)\,=^=-Dj{-x)U=a-l. (4.17) 
s>n ^ ' 

By the definition of D it follows that -D[/(— a;)] |^^^ = —D\f{x)\\x=-c-\- Repeated 
application of this implies 

^"[/(-^)]U=c = (-l)"i^"[/(a;)]U=_e-n. 
Using this, the right hand side of (4.17) becomes 

(-1)' 



On the other hand, the left hand side of (4.17) equals 



Therefore, (4.17) is equivalent to 

(~''r'')i?"+7(^)U=l = i?'^/(x)U=i_e-n. (4.18) 

i>o V / 

However, this is readily deduced from Newton's formula for a polynomial j{x) in terms 
of powers of the difference operator (see e.g. [13, (4), p. 75]): 

/(^) = /(rf) + ^ '^)^/(^)U=d + ("^ 2 + • • • 

Indeed, taking x — \ — c — d — \ m. the above equality, and applying then to both 
sides, we obtain equality (4.18). This completes the proof of (4.12). □ 

We are now ready to present the proof of the main result of this section. 
Proof of Proposition 4-1- By [14, (4.11)], for a; < — 1 the function P given by (1.3) can 
be written as 

P{x,y) = -^ I t-y-\-l-t)—'dt, (4.19) 
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where the integral is a contour integral along the the counterclockwise oriented arc of the 
unit circle connecting e^'^*/^ to e^'^*/^. Therefore, for all i? > a we have 



^4,ri/3 

^ I -bt 1 4-\R-a. 



P{-R - 1 + a, -1 + 6) = -— / ^-"(-1 - t)''-''dt 



27T 

i 

'2^ 



g27ri/3 
4Tvi /3 



27ri/3 




/3 



Je47ri/3 



(4.20) 



where p{t) and (5(t) are given by (4.2)-(4.3). We find the asymptotic series for large R 
of the two integrals in the curly braces using Laplace's method for contour integrals as 
presented in [19, Theorem 6.1, p. 125]. 
Denote the two integrals by 



Ii{R)= f e-^P^'^Q{t)dt 

7g27ri/3 

I^{R) = / ' e-''P^'^Q{t)dt. 



One readily verifies that Ii{R) satisfies conditions {i)~{v) of [19, pp. 121-122]. Indeed, 
p{t) and Q{t) are clearly independent of R and holomorphic in the open disk D(e^'^*/^, 1), 
which contains the interior of the integration path V; this checks conditions (i) and (ii). 
Condition (iii) requires p{t) and Q{t) to admit power series expansions around t = e^^*/^; 
this is clear from their definitions. Condition (iv) requires Ii{R) to converge at —1 
absolutely and uniformly for R large enough; this is readily seen to be the case. Finally, 
to check condition (v) we need to show that Re{p{t) -p(e^'^*/'^)} > for t e (e^'^*/^, -1). 
Since Rep(t) = Re{— ln(— 1 — t)} = — In | — 1 — t|, and | — 1 — t| < 1 throughout V, with 
equality only for t = e^'^*/^, the latter requirement is verified as well. 

Therefore, [19, Theorem 6.1, p. 125] applies and the asymptotics series of Ii{R) is given 

by 

S=0 \ r / 
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where ( = e'^'^'^^^, jj. is the smaUest power of t — z with a non-zero coefficient in the series 
expansion of p{t) —p{C)j A — 1 is the smaUest power of t — with a non-zero coefficient in 
the expansion of Q{t), and is the coefficient of in the power series expansion (4.12) 
of Q{t)/p'{t). Since the value of po in Lemma 4.3 is p'{C) = —1/(1 + C) 7^ 0, we have 
11=1. Lemma 4.4 shows that A = 1, and also provides the explicit form of the a^'s. 
Therefore (4.21) becomes 



p--Rp(C) 



s=0 



E 

.1=0 



=o-|-6-lC" 



(4.22) 



The sum in the brackets is, by the binomial theorem, just (1 — D(^ ^) ^{x'^)\x=a+b-i- 
Since -1 - ( = (-^ and e"'^^'^^) = (-1 - C)^ = C"'^, we obtain from (4.22) the explicit 
asymptotic series of Ii (R) to be 



h{R) ^a-b-l ^ (l-i^ C-^)-^(^")U=a+6-l 
s=0 



c- 



(4.23) 



The value of the integral l2{R) is readily seen to be the complex conjugate of Ii{R). 
It follows then from (4.23) that 



h{R) +5+1 V a - DO-H^')\.=a+b-l 

^ s=0 



(4.24) 



Substituting (4.23)-(4.24) with R replaced by 3r in (4.20) and using = 1 we obtain 
(4.1). □ 



5. Asymptotic singularity and Newton's divided difference operator 



Direct application of Proposition 3.2 to the collection of + Si=i triangular 

holes of side 2 that the multiholes E^^ , . . . , E^^ , W^^ consist of leads to the 

following result. 

Let S = Xlili Si and T = tj. 
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Lemma 5.1. For S >T, the correlation on the left hand side of (2.1) can be written as 
Lb{El^{Rxu Ryi), E^jRxm, Rym), W^^{Rzi, Rwi), W^^iRzn, RWn)) 

= |det(M)|, 

where 

M^[C^ ■ ■ ■ Cn D] (5.1) 
is a 2S X 2S matrix whose blocks are given by 



P(-R(zj-xi)+ aii-bji-l, -R(wj-yi)+ g(an- bji)-l) 



P(-ie(Zj-J:,)+a,,- 



-R(^j-yi)+<l(ai\- 



P(-R(zj- xi)+ Oj J - bji - 1, -R(wj-yi)+ qiPi^-bjO-l) 



P(-R(zj-Xi)+ai,-bj,-l,-R(irj-yi)+q(aurl'j,,)-l) 









P{-R(.Zj- ^„ ) + 1- ijl - 1 , ym)+ 1- *Jl ) - 1) 






P(-R(zj- x„)+a„i- bj,.- 1 , -R(wj-y„)+ q(a„ i,;, ,)- 1) 


















P(-R(zj-x„) + a„s- bji-l, -R{wj-y„)+ q(.am,-bji)-l) 






P(-R(zj-x„)+a„,-bj,.-l,-R(wj-yJ+q(a„s„-bj,^)-l) 



















(5.2) 



for j = 1, . . . ,n, and 
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Uo(Rxi+ ail, ^yi'^ 1^11^ 









Vs_j_l (Rxi+ flij, Ryi+ qaii) 



















D = 



UoiRXm+Oms^' Rym+^amsJ 






f/s_7--l(^^m+«m*™' Rym+^amsJ 















(5.3) 

(Cj has 2S rows and 2tj columns, D has 2S rows and 2S — 2T columns); here for a 
function f of two arguments the box 



f(x,y) 









denotes the 2x2 matrix 



f{x,y) f{x-l,y+l) 
f{x + l,y-l) f{x,y) 



and the large boxes with rounded corners outline 2sj x 2tj blocks of (3.4) that correspond 
to a given pair {E^ -, WiyJ of multiholes. 
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Therefore, in order to determine the asymptotics of the correlation of the multiholes, 
it suffices to find the asymptotics of det(M) as ^ oo. 

In working out the asymptotics of the determinant of a matrix whose entries are 
functions of a parameter two situations can occur. Consider the matrix formed by the 
asymptotic approximations of the entries of the original matrix. If this matrix is non- 
singular, then its determinant gives the asymptotics of the determinant of the original 
matrix. On the other hand, the matrix of approximants may be singular. Then to find the 
asymptotics of the determinant of the original matrix, one needs in general to work with 
higher order terms of the asymptotic series of the entries, and both finding these higher 
terms explicitly and evaluating the resulting determinant gets usually very complicated. 

It turns out that our matrix M falls into the second category: the matrix of the 
asymptotic approximants (as R oo) of its entries is singular unless si = • • • = Sm = 
ti = ■ ■ ■ = tn = 1- 

Fortunately, however, it also turns out that we can find suitable determinant-preserving 
row and column operations on M that transform it into a matrix M' which falls into the 
first category. The row and column operations that achieve this are patterned after the 
operator involved in the following result. 

Lemma 5.2. Let X he an n x n matrix, and suppose rows ii, . . . ,ik in it are the vectors 
/(ci), .... f{ck), respectively, where f is some vector function. Consider the operator fx 
that acts on rows zi, . . . , of X by transforming them as 



r/(ci)i 




r ^^V(ci) 1 


/(C2) 


1— > 


(C2 - Ci)PV(ci) 






-(Cfe -Ci){Ck -C2)...{Ck - Cfe_i)P^-V(ci). 



where V is Newton's divided difference operator, whose powers are defined inductively 
by V^f = f and V^fia) = {V^-^f{ci+t) - V^'^ f (a)) / {a+r - Ci). Then det(/x(X)) = 
det(X). 

Proof. By [13, (2), p. 20], for any function / and any sequence ci, . . . , q one has 

/(q) = /(ci) + (q - ci)P/(ci) + (q - ci)(q - C2)V^f{ci) + ■■■ 

+ {Cl - Ci)(q - C2) • ■ • (q - Q_i)P'-V(ci). 

(5.5) 

Define a sequence of matrices Xj as follows. Let Xq = ^(X); let Xi be the matrix 
obtained from Xq by replacing row ^2 by the sum of rows ^2 and ii; let X2 be the matrix 
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obtained from Xi by replacing row by the sum of rows Z3 and i2', and so on. Clearly, 
rows ii, . . . , ifc of Xk-i are then i?^^, Ri^ + Ri^, . . . , Ri^ + ■ ■ ■ + Ri^, respectively, where 
Ri stands for row i of //(X). By (5.5), these are exactly rows ii,...,ik of X. Thus, 
= X. Since X^-i was obtained by a sequence of elementary row operations from 
//(X), it follows that det(//(X)) = det(Xfc_i) = det(X). □ 

We now define the matrix M' that resolves the asymptotic singularity of M described 
before the statement of Lemma 5.2. 

By its definition in the statement of Lemma 5.1, it is apparent that for any i = 
1, . . . , m, rows 2 + 1, 2 Sjt + 3, . . . , 2 X1a;=\ -^fc + 2si - 1 of M are of the form 

f{0'ii): f{0'i2), ■ ■ ■ , fi^isi), for a Suitable vector function /. Apply the operator (5.4) on 
these Si rows of M, for each i = 1, . . . ,m. 

By the same token, rows 2 Ylk^i Sk + 2,2 Yj^i—i ■Sfe + 4, . . . , 2 Sk + 2si of M are 

also of the form /(aji), f{o-i2)i ■ ■ • ifif^isi)-, for suitable vector function /. Apply the 
operator (5.4) on these Sj rows of M as well, for each i = l,...,m. Let Mi be the 
resulting matrix. 

In the same fashion, for any fixed j = 1, . . . , n, it is clear from (5.1) that the columns of 
M with indices 2 Yli'i=i + 1? 2 ^^1=1 + 3, ■ • • , 2 Yli'i=i + — 1, on the one hand, and 
those with indices 2 Yl'lll ti + 2,2 Yl'IZi + 4, . . . , 2 Yl'lll ti + 2tj, on the other, are each 
of the form f{—bji), f{—hj2)i . . . , f{—bjtj), for a suitable vector function /. Furthermore, 
the same is true for the corresponding two tj-tuples of columns in matrix Mi. We define 
M' to be the matrix obtained from Mi by applying the operator (5.4) — which clearly 
works equally well for columns — along each of these tj-tuples of columns, for j = 1, . . . , n. 
Thus M' is obtained from M by applying operator (5.4) for a total of 2m + 2n times. 

To describe the entries of M' the following terminology will be convenient. We refer to 
the first 2T columns of M as the P-part of M; the last 2S — 2T columns form the U-part 
of M (see (5.1)-(5.3)). The P-part consists of an m x n array of blocks, the block Pij in 
position {i, j) being in turn a x tj array of 2 x 2 matrices^. The t/-part of M is a column 
of m blocks, the block Ui in position i being an Si x {S — T) array of 2 x 2 matrices. 
Block-row i consists of Pn, Pi2, ■ ■ ■ , Pim, Uf, block-column j consists of Pij, P2j, . . . , Pnj 
(it is precisely this block-column that is displayed in (5.2)). The 2x2 block-matrices 
naturally group the 2tj columns of block-column j into tj pairs of consecutive columns; 
we call them bi-columns; bi-column /c of a block-column consists of columns 2k — 1 and 
2k of that block-column. Bi-rows are defined analogously. 

Since the block structure of Mi and M' is clearly the same as that of M, we use the 
above terminology for them as well. 

^The block can be thought of as corresponding to the pair (-El - , W? ) of multiholes; then its 2 x 2 
block in position {k,l) corresponds to the pair {E{aik,qaik),W{bji,qbji)) of triangular holes of side 2. 
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With this terminology, a concise description of our constructions is that we first obtain 
Ml from M by applying operator (5.4) along the top halves of the bi-rows in block-row 
i, and then along the bottom halves of the same bi-rows, for all i = 1, . . . ,m. M' is 
obtained from Mi by applying operator (5.4) along the left halves of its bi-columns in 
block-column j, and then along the right halves of the same bi-columns, for j = 1, . . . , n. 

The 2x2 matrix in position {k, I) in the (z, j)-block of the P-part of M (1 < A; < Sj, 
1<1< tj) is 



P(-R(zj - Xi) + ai;.-bji-l,- RXw j -yi) + g(aii.-bji)-l) P{- R{zj — Xi)+aii^-bji-2,- R(-wj - J/i )+9(aifc - bj( ) ) 

Pi-R(zj -x^) + a^f.-bjl,-R('Wj-v^) + q(a^f.-bjl)-2) P(-R(zj - x^)+a^f.-b ,- R^w j - y^)+q{,a^f.-bjl)-l) 

By our construction, the corresponding 2x2 matrix in M' is 

-P( — -R(zj-a:j) + aji-()ji-l, — -R(«)j— j/i) + gr(aji-()ji)-l) P( — R{zj - x j ) + a ji — i — 2, - - yi)+i3(aji -bji)) 

P(-R(zj-x^) + a^■^^-bj'^^,-R(■Wj-y^) + q(a^-^^-bj-i^)-2) P(-R(zj-Xi) + aii-bji-l,-R(wj-yi)+g(a^i-bji)-l) 



(5.6) 



where the powers of T> act entry- wise, ^ acting on the sequence an, aj2, 
on the sequence —bji, —bj2, ■ ■ ■ , —bjtj , and 



l-l 



a 



k 



{aik — aii){aik — 0^2) • • • {aik — ai^k-i), 



Pi = i-bjl + bji){-bji + bj2) ■ ■ ■ {-bji + b^,i-i). 



(5.7) 
(5.8) 



On the other hand, the 2x2 submatrix of M at the intersection of rows 2k — 1 and 
2k of block-row i with columns 2r + 2Z - 1 and 2T + 2Z (1 < /c < s,-, 1 < Z < 5 - T) is 
by (5.1)-(5.3) 

Ui-i {Rx^ + ttifc, Ryi + qaik) Ui-i {Rx^ + a^fc - 1, Ryi + qaik + 1) 
Ui-i{Rxi + aik + 1, Ryi + Q^ik - 1) Ui-i{Rxi + a^k, Ryi + qa^k) 

By our construction, the corresponding 2x2 submatrix of Mi — and thus of M', as our 
applications of the operator (5.4) to Mi do not affect the t/-part of Mi — is 



Ui-i{Rxi + an, Ryi + qan) Ui-i{Rxi + an - l,Ryi + qan + 1) 
Ui^i{Rxi + an + 1, Ryi + qan - 1) Ui-i{Rxi + an, Ryi + qan) 
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(5.9) 



with a)^ given by (5.7). 

Thus in short, the 2x2 blocks that form M' — the matrix of the same determinant 
as M, but whose entries' approximants wiU turn out to form a non-singular matrix — are 
given by (5.6)-(5.9). 

The asymptotics of the entries of (5.6) and (5.9) is worked out in the next two sections. 
The formulas we obtain there lead us to the following result. 

Proposition 5.3. The asymptotics of the determinant of the matrix M in the statement 
of Lemma 5. 1 is given by 

\2S m n 

^ n n K-«ifc)'n n (K-bik? 

^ i=l \<i<k<ai i=l l<j<k<ti 

X det(iW^") i?^'f2I)i<i<j<m ^■i*j+5Zi<i<j<T!, 5Z"=i Si^j} 

(5.'l0) 

where 

M" = [C" D"] (5.11) 
is the 2S x 25' matrix whose blocks are given by 



a"(0,0,zi-a:i ,\rj-yi ) 






a"(0,(i - l,zi-*i ,ivi-> i ) 


a'2 


• • • 


a"(0,0,z„-i1 .Wn-yi ) 


a'^ 




a"(0,(„- l,Z„-il .Wb-JI ) 


a'2 




a" 


A^' 


a" 


A" 


a" 


A^' 


a" 












} 


a"(ii -l,0,Zi-:ti ,w,-J| ) 


a'2 




a"(ii -1,(i-1,Zi-a:, ,wi-) | ) 


a'' 


a"(si -1,0,Z„-;ci .Wn-J'i ) 


a'2 




a"{.Si - l,Zn-J:i ,w„->i) 


a'2 


a" 


a" 


A^' 


a" 


A^l 


a" 


A^' 


a" 


• 
• 
• 




• 
• 
• 


A"{0,0,z,-x,n.w,-y„) 


a'2 




a"(0,(i-1,zi-i„,wi-j„,) 


a'^ 


• • • 


a"(0,0,z„-j„,w„-j,„) 


a'2 




A"lO,tn-l,Z„-x„,Wi,-y,n) 


A>^ 




a" 


A^> 


a" 


A^l 


a" 


A^' 


a" 














A^\s„-l,0,Zi-x„,Wj-y„,> 


a'^ 




a"(j„-1,(i-1,Zi-a:„,wi-}',„) 


a'^ 


a"(s„-I.O,z„-j,„.w„-j,„) 


a'2 




a"(i,„ - l,t„- l,z„-J:,„,w„->i„) 


a'^ 




a" 




a" 


A^> 


a" 


a" 


a" 



and 
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(5.12) 



b"(0,*i,>.,) 




*0 







"s-T-l 


B^' 

"S-T-l 


b" 

"S-T-l 



fii'ta - 1, 


"a 






r'^ 


r" 


b" 




B^' 


"s-r-1 


• 
• 
• 


«"(0,j„,.y,„) 


p 

Bo 




B]ij._,(0..i,„,v,„) 


r'2 

^5- T- 1 


b^' 






B^-r-i 


r" 

"S-T-l 












Bo 




Bjij-.i (J„ - 1, Im.J,,,) 


r'2 




Bq 




Bj-r-i 


„" 

B^-r-i 



(5.13) 



TO 



th^ 



11 
12 
21 



k + l 
k 



-1 (i-gC)^+' 



(i-gC-^)fe+' 

(5.14) 



and 



5, 



11 



12 
21 



{k,u,v) 



-l\l-k 



(5.15) 



^Equality (5.14) defines simultaneously three functions, A^^, A^^ and A'^^; for instance, to obtain 
the expression for A^^, one chooses the middle quantities in the curly braces on the right hand side. 
Similarly in (5.15). 

32 



{here ( = e^'^*/^; in (5.12) and (5.13) the arguments within each 2x2 box are the same; 
for brevity, they are indicated only for the top left entry) . 

Proof. By Proposition 7.1, as i? — > oo, the entries of the 2x2 matrix (5.6) have their 
asymptotics given by 
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A^^ (k — l,l — l,Zj —Xi ,Wj —yi) A^'^ (k — l,l — l,Zj —Xi ,Wj—yi ) 



A^'^{k-l,l-l 



-yi) A^\k-l,l-l 



-Vi) 



+ 0{R-^-^). 



Therefore, the asymptotic approximations (as R — > oo) of the entries in the (^, j) block 
of the P-part of M' form the matrix^ 



VET. 

27r 



■X 



a^'^^P^ All (0,0,2,- _^ . _y,^ 


(i 




a^*^/3t'^'A"(0,tj-l,2i-x,-,Wi-yj) 




R 




R • • 

^0)^12 






R 




R • • 






ai'^l3['^A^\l,0,Zi-Xj,Wi-yj) 


(i 
0=2 


(3[''^A^^ 








(i 


R^ • • 






fl2 




R2 • • 




R'J+' 



,(p[JU^\s,-l,0,z,-xj,Wi-y,) gffM^'^^' 



R^'i 



(5.16) 

and the difference between each entry of M' and its approximant shown in (5.16) is 
0{R^~^), if the order of the approximant is R'^ (as in (5.12), the arguments of A^^, A^^ 



^To distinguish it from a variable index, we denote in this section the complex number i by y/—! 
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and A"^^ arc the same within each 2x2 block in (5.16); to economize space, they are 
indicated only for the top left entry in each such block). 

On the other hand, by Proposition 6.1, matrix (5.9) can be written as 



a 



27r R''-^ 



Thus, the asymptotic approximations of the entries in the z-th block of the t/-part of 
M' form the matrix 



27r 



R^i-^ 



R" 



-(S-T) 
,21 

'S-T-1 











-Ds_T-l 


RO 


i?o • • 






i^l-(S-T) 


"1 -"0 








ai 'Bs_T-i 


RO 


RO ■ ■ 


JJl-(S-T) 






ai'^Bl\l,Xi,yi) 






,yi) 


_ (i) r)12 

"^2 -°S-T-1 


R^ 




fl2-(S-T) 




^2-(S-T) 










«2 -Bs-T-l 


Ri 


Jti • • 


^2-(S-T) 




^2-(S-T) 



„ (<) d12 



-(S-T) 
'S-T-1 



R^ 



-(S-T) J 



(5.17) 

with the same conventions about arguments and the same statement about the difference 
between exact values of entries of M' and their above approximations as in (5.16). 

The forms (5.16) and (5.17) of the blocks of M' clearly show that all the factors 
— V— l/(27r) factor out along the rows, for a total contribution of (l/27r)^'^ to det(M'). 
It is also apparent from (5.16) and (5.17) that all the cu's factor out along the rows of 
M', giving, by (5.7), a combined contribution of 



Tlx S ' TTh 



i=l fc=l 



1=1 l<i<k<Si 
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In a similar fashion, all the /?'s factor out along the columns of M', and by (5.8) they 
combine to a multiplicative factor of 

n tj ^ n 

=n n c-^.-M^ 

J = 1Z=1 1=1 l<j<k<ti 

Furthermore, the factors in R can also be factored out along the rows and columns of 
M'. Indeed, multiply the rows in block-row i by i?^, i?^, i?^, i?^, . . . , respectively 
(from top to bottom), for z = 1, . . . , m. By (5.16)-(5.17), this will make the degrees in R 
of the entries constant along the columns. More precisely, in block-column j, the degree 
in R becomes along columns 1 and 2, —1 along columns 3 and 4, and so on, reaching 
degree — (tj — 1) along columns 2tj — 1 and 2tj of block-column j; and in the [/-part the 
degrees in R become 1 for the first two columns, 2 for the next two, and so on, reaching 
degree 5 — T for the last two columns of the [/-part. Factoring out these powers of R 
along the columns, and taking into account the factors wc multiplied the rows by, one 
obtains that the overall exponent of R that gets factored out of the determinant of M' is 

m Si n tj S — T m n 

i=l k=l 3=1 1=1 3=1 i=l 3 = 1 

which is readily seen to equal 2 YliiKKjKm ^^^j <i<j<n ^i^3 2 X^i=l X^j=l ^i^j- 

However, after factoring out this way the — V— l/(27r)'s, cu's, /?'s and the powers of R 
from the blocks (5.16) and (5.17) of M', the resulting matrix is precisely the matrix M" 
given by (5.11)-(5.15). This completes the proof. □ 

6. The asymptotics of the entries in the U-part of M' 

In this section we determine the asymptotics of the entries of (5.9) as i? — > oo. These 
asymptotics follow from the following result. When Newton's divided difference operator 
V acts on a function / of more than one variable, we write Vxf to indicate that V acts 
on / regarded as a function of the variable x. 

Proposition 6.1. Let q e Q satisfying 3\1 - q. Then for any fixed c,d e Z, u,v e 3Z, 
and < A;, Z e Z, we have 

V^Ui {Ru + x + c, Rv + qx + d) |a;=ai = 

- ^ Q [^-'^-'(1 - qc,)\u - vcy-^ - (-^=+^+1(1 - qc^f{u - vcy-"] R'-" 

+ 0(i?^-^-i), (6.1) 
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where C = e^'^*/^ and acts with respect to a fixed integer sequence ai, 02, . . . with the 
property qaj e Z, j > 1. 

In our proof of the above result we use the following three preliminary lemmas. 

Lemma 6.2. For any a,b & 'Z and < k & Z, we have 

(l—D(^~^)~^(x^)\x=a+b-i — io,—b()^+ monomials in a and b of joint degree < k, (6.2) 

and 

{l—D()~^(x'^)\x=a+b-i — {a'—bC~^)^+ monomials in a and b of joint degree < k. (6.3) 

Proof. By repeated application of the fact that {x + 1)'="'^ — x''~'^ = {k — n)x''~'^~^ + 
as a; — > 00, one obtains that 

D^'x'' = k{k-l)---{k-n + l)x''-'' + 0{x''-''-^), x^oo. (6.4) 

Since 

(1 - Dc'r'ix'^) = x^- [\y~'Dx^ + ("2^)^-'^'^' - • • • , 

we get by (6.4) that 

(1 - DO-\x'^)\,=a+b-i ={a + bf + bC^k{a + bf-' + ^C-'A;(A; - l)(a + bf'^ 

+ %C''k{k-l){k-2){a + bf-^ + 

+ monomials in a and b of joint degree < k (6-5) 

(we used here that the portion of {^^'^ of maximal degree is i,-). However, the leading 
part on the right hand side in (6.5) can be written as 

(a + hf + (J) (a + br-'bC' + Q (a + bf-\bC'f + ■ ■ 

= [{a + b)+bC-'f 
= {a-bO\ 

since 1 + C,~^ — Thus, (6.2) follows from (6.5). Replacing by ( in the above 
argument proves (6.3). □ 
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Lemma 6.3. Let hn{xi, . . . , Xs) be the complete homogeneous symmetric function of de- 
gree n in the variables xi,...jXs {hn := for n < 0, ho = 1, hi = Ylk^k, ^2 = 
Tlik -^k '^k<i ^kXi, and so on; see e.g. [25, p.294]). Then for any < n, k & Z we have 



1> X''\x=ai = hn-k{ai, . . . , ttk+l). (6.6) 

Proof. We proceed by induction on k. For /c = 0, P x^ I x=ai — '^1^ ) which, is cle&rly the 
same as hn{ai). The induction step follows by the calculation 



Ti^^'^l — \x=a2 ^ \x=ai 

X \x=a\ — 



^n-A;+l(02, • • • , Ofc+l) — ^n-A;+l('^l) • • • ) ^^k) 



n— fe+1 T 



= E 



fe+l^n-fc+l-r(Q^25 • • • 5 (ife) — 0ii/t„_fe+i_r(a2, • • • , Ofe) 



r.=0 '^'=+1 ~ ^1 



n—k 



hn-k-r{0'2, • • • , Ofe)(Oi + ^^fe+l + • • • + a^_,_i) 



r=0 



— hn-kio-l-, • • • ) Olfe+i), 

where at the second equality we used the induction hypothesis. □ 
Lemma 6.4. For any <k,l eZ and A, B,C e C we have 

Vl{Ax + RB + C)%=a, = {^^A^B^-^R^-^ + 0{R^-^-'^), R^oo. (6.7) 

Proof. Expand the power on the left hand side by the binomial theorem and use (6.6) 
to obtain 

Vl{Ax + RB + Cy U=a, = 5^ ) A-(i?S + Cy-^hn-k{ai, ak+i) 

\k T-fl — k -rtl — k u /„ „ ^ 1 /O/D^ — fc— 1> 



n=0 
'I 
k 



A'B'-^R'-^ho{ai, ak+i) + 0{R'-^-'). 



Since /iq = 1, this completes the proof. □ 
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Proof of Proposition 6.1. By (4.1) and Lemma 6.2 we obtain 

Ui{a, h) = - ^[C-'-\a - bCY - C-"+'+H« - ^'C-')^] 

+ monomials in a and b of joint degree < I. (6-8) 

In particular, 

Ui{Ru + a; + c, Rv -\- qx -\- d) 

_ ^_R^u-v)-x{l-q)-c+d+l^(^-^ _ ^^-1)^ ^ _ ^^-1) ^ ^ _ 

+ Yl Ca,f3{R'^ + ^ + cr{Rv + qx + d)^, (6.9) 

a,;3>0 
a + 0<l 

where the Ca^p^s are independent of R and x. 

By hypothesis, R{u — v) + x{l — g) is a multiple of 3. Thus, since C = e^'^*/^, the 
exponents of C in front of the square brackets in (6.9) simplify to c — d — 1 and —c + d + 1, 
respectively. 

By Lemma 6.4, as — > oo we have 



P^[(l - qC)x + R{u -vQ + c- dCY\a.=a, = Q (1 - qOH^ - vO'-'^R'-^ + 0{R'-'^-^ 



) 



and 



^'[(1 - ^r')^ + R{u - vc') + c - dc']%=a, = (!) (1 - 9C-')'(« - vc^y-^R'-^ 



+ 0{R^-^-^). 



On the other hand, the result of applying to the sum in (6.9) is clearly 0(i?' ^ ^). 
Applying to both sides of (6.9) and using the above observations we obtain (6.1). □ 

7. The asymptotics of the entries in the P-part of M' 

In this section we determine the asymptotics of the entries of (5.6) as — > oo. The 
resulting formulas are contained in the following result. As in the previous section, a 
variable at the index of Newton's divided difference operator V indicates that the function 
upon which it acts is regarded as a function of that variable. 
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Proposition 7.1. Let q e Q satisfying 3|1 - q. Then for any fixed c,d e Z, u,v e 3Z, 
(u, v) ^ (0, 0), and < k,l & Z, we have 



V[ {DlP{~Ru + x + y + c, -Rv + q{x + y) + d)\,=a, } 

i fk + l\ 
~2^V I J 



y=bi 

l\k+l 



" (7.1) 



where ( = e^'^^l'^ , acts with respect to a fixed integer sequence ai, 02, . . . , T>y acts with 
respect to another integer sequence 61, 62, ... , and qaj, qbj e Z for all j > 1- 

In our proof we will employ the following preliminary results. 

Lemma 7.2. For < n,k & Z and indeterminates q, z, cq, ci, . . . , one has 



k {-iy{qco)j (co)fc-j ^k-2i {-iy{qci)j {ci)k-j ,fc 



EK 
i=o 



71 {k-j)r ^--3 



{k-jy: 



-2i 



(co - ci)(co - C2) ■ • • (co - Cn) (ci - co)(ci - C2) • ■ • (ci - c„) 



+ 



+ 



Y-fc (-l)^(9C„)j (c^)fc-j ^fc_2,- 
2.,=o JT (fc-j)!^ 



0, if A; < n, 

(Cn - Co)(c„ - Ci) • • • (C„ - Cn-l) ] i.^^ _ g^"!)"^ if /c = n. 



Proof. Multiplying by no<j</<n('^' ~ ^j)^ statement becomes equivalent to the 
polynomial identities 

(-1)" n ('=.-.)E^=^4r^|^-'-^' + - 



0<j<l<n 



J=0 



+(-i)° n («--.)E^=%^7F^-'-^^ 

0<j<l<n j=0 •' ' ^ 

0, if A; < n, 

= < n (^'-^^•) n . . 

0<.</<n .Y^(^ Y_iyq3^n-2j^ H k = Tl. 

3=0 ^-^^ 
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(7.2) 



(7.3) 



We prove (7.3) first, so let k = n. Agreement of the coefficients of z'^ on the two 
sides of (7.3) amounts to 



(-1) YY (q-C,) — ^-^ + ... 



s\ in — s)\ 

0<j<l<n ^ ' 



+ (-1) 11 (Q-c,) 



{n-s)\ 

0<j<l<n ^ ' 



1 / n 



Q(-i)v n (^'-^i)- 



0<j<l<n 

Simphfying this identity we obtain that (7.3) is equivalent to the set of identities 

(-l)"(gCo)s(co)n-s H {ci-Cj)-\ \-{-lf{qCn)s{Cn)n-s H {ci - Cj) 

0<j<l<n 0<j<l<n 

=Q' n (^'-^^•)' (7-4) 

0<j<l<n 

for < s < n. 

Let < m < p < n, and assume = Cp. We claim that this causes the left hand side 
of (7.4) to be 0. Indeed, let A(xq, . . . ,Xt) = Y[o<j<l<ti'^i ~ "^j)' Then for = Cp the 
left hand side of (7.4) becomes 

( 1) ^(cQj • • • ) C^n—lj Crffi^i, . . . , Cp, . . . , C^) (?C^)s (c^),^_s + 

( 1)^^(C0; • • • ) CjTj, . . . , Cp—i, Cp-|-i, . . . , Cn^(^QCp^g(^Cp^ri—s- 

However, since Cm — Cp, the second list of indeterminates on which A acts above is 
obtained from the first such list by moving Cp across p — m — 1 neighbors to its left. This 
implies that the two A-expressions differ by a multiplicative constant of (— 1)^""*"^, and 
thus the above expression equals 0. 

Regard the two sides of (7.4) as polynomials is Z[q'][co, . . . , c^]. As a consequence of 
our observation in the previous paragraph, the polynomial on the left is divisible by the 
product no<j<z<n('^' ~ ^d)- Furthermore, the degree of the polynomial on the left is 
(2) ^' "which is the same as the degree of this product. It follows that the two sides of 
(7.4) are equal up to a multiplicative factor in 
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To finish proving (7.4) it suffices to show that it holds for a speciahzation of the c/s. 
For Cj = j, j = 0, . . . ,n, this amounts to 

J2 f -M ., id)s{j)n-s = q', 0<s<n. (7.5) 

Take f{x) = and w = in the expression (4.16) for the n-th power of the divided 
difference operator D. We obtain 

( 1 f 0, if < t < n, 

Thus, only the term in in the expansion of 

{(i3)s {j)n-s = qj{qj + 1) • • • (^j + s - l)j(j + 1) • • • (j + n - s - 1) 
= q^j'^ + terms in j of degree < n 

contributes to the left hand side of (7.5). This proves (7.5), and completes the proof of 
(7.3). 

We now turn to proving (7.2), so assume k < n. By extracting the coefficients of 
^fe-2s^ (7.2) becomes equivalent to 

(-l)"(gco).(co)fe-. n {ci-Cj) + ---{-lf{qCn)s{cn)k-s n = (7-7) 

0<j<l<n 0<j<l<0 

for < s < A;. 

By the arguments we used in the case k = n, the left hand side of (7.7) is seen to be of 
the form a Y[o<j<i<n(^i ~ ^3)1 with a G Z[g]. To establish (7.7) it suffices to show ct = 0. 
The specialization that led to (7.5) shows that in order to deduce the latter it is enough 
to check that ^ 

i~^K!-M j)^ (j)k-s = 0, < s < A; < n. 

This follows from (7.6) by the same reasoning we used to prove (7.5). □ 

Lemma 7.3. Let q E Q with 3|1 — g and consider a sequence ao,ai, . . . of integers so that 
qaj e Z, J > 0. Then for any n> 0, in the power series expansion 

V^t-^^{-l-t)-%=ao=J2^k{t-0' (7.8) 

fc>0 
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around C, = e^'^*/^, the first non-zero coefficient is Qn = ^{C ~ qC ^)^- 
Proof. By [13, (l),p.l9], for any functfon / we have 

V^f(ao) = + + 

(ao - ai)(ao - 02) • • • (oq - a^) (ai - ao)(ai - 02) • • • (ai - a^) 



(an - ao)(an - ai) • • • (an - an-i) " 

(7.9) 

Take /(x,t) = t~^^(— 1 — t)~^ in the above equahty. Then the resulting expression wiU 
supply the power series expansion (7.8), provided we work out the coefficients of the series 
expansion 

fc>0 

Let g{t) — t^*^ and h{t) = (— 1 — t)~^. Their derivatives of order / are given by 

g^'\t) = {~iy{qx)it-'i--', 
h^i){t) = {x)i{-l-t)—'. 
Thus, by Leibniz's formula for the derivative of a product. 



= E (*)(-i)'to%),K).-.C-«-'(-i - 0- 

1=0 ^ ^ 



■aj-{k-l) 



1=0 

where at the last equality we used — 1 — C = (^~^ and (^"^(^"9) = i Substituting this in 
(7.9) we obtain that the coefficient of {t—Q^ in the power series expansion of V^t~'^^{—1 — 
t)~^|x=ao around t = ^ is 

v-fe (-l)^(Qao)j {ao)k-j 



(ao - ai)(ao - a2) • • • (ao - an) 



(-l)^(gai)j {ai)k-j .k-2'i v-fe {-iy{qan)j Mk-j .k-2'i 



(ai - ao)(ai - 02) ■ ■ ■ (ai - a^) (a^ - ao)(an - ■ ■ • (a^ - a^-i) ' 

Apply Lemma 7.2 to complete the proof. □ 
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Lemma 7.4. Let g e Q with 3|1 — g, and let ai, 02, . . . and bi, 62, . . . be sequences of 
integers so that qaj, qbj & Z, j > 1. Then for any k,l >0, in the power series expansion 

vl{v',t-'^^^+y\-i-t)-^^+y^U=aA\ ^ =$]gs(t-cr 

around ^ = e^'^*/^, the first non-zero coefficient is qk+i = -^(C ~ ?C~^)'°^'- 
Proof. We have 

I, J l?/=6i 



|X=Ol • 



The statement foUows therefore by Lemma 7.3. □ 
Proof of Proposition 7.1. Assume first that u> 0. Then by (4.19) we obtain that 

P{-Ru + x + y^-c,- Rv + q{x + y) + d) 

27r Jg27ri/3 

1/ ' e-^P^'^Q{t)dt- f ' e-^PWg(t)d4 



I 

'2^ 



for all sufficiently large R, where 



p{t) = -v\nt-u\n{-l -t), (7.10) 



Applying T>!yT>^ to both sides yields 



Vl {VlP{-Ru + x + y + c,-Rv + q{x + y) + d)U=„, } |^^^^ 

I / e-^P^'^Q{t)dt - I e-^P^'^Q{t)dt\ , 



'2^ 



(7.11) 
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where p{t) is given by (7.10) and 

Q(t) = t-'^-i(-l-t)-^-iDj^ |D^t-9(-+?^)(_l_t)-(-+^)|^=„ 1 . (7.12) 

As in (4.20), the values of the two integrals in (7.11) are complex conjugates. Let 



IiR)= [ e-^P^'^Q{t)dt. 



We find the asymptotics of I{R) by the approach we used in the proof of Proposition 
4.1 — Laplace's method for contour integrals as presented in [19, Theorem 6.1, p. 125]. The 
integral I{R) is readily seen to satisfy conditions {i)-{v) of [19, pp. 121-122] {{v) holds 
since we are in the case u > 0). 

The first two coefficients of the power series expansion 

p{t) = p{C) +Po{t-C)+Pi{t-C? + --- 

of the function p{t) of (7.10) are easily checked to satisfy e'^P^'^^ = = 1 (the 

second equality since u and v are multiples of 3) and po = p'{C) = u( — Since 
Po 7^ it also follows that the value of fx when applying (4.21) for I{R) is /i = 1. 

The function Q{t) given by (7.12) is precisely the function whose power series expansion 
is considered in Lemma 7.4, multiplied by t~'^~^{—l — t)~^~^. Since the first non-zero 
coefficient in the power series expansion of t^"^^^ {—1 — t)~'^~^ around t — ^ is the constant 
term, and it equals — C)~^~^ = C^~^i we obtain by Lemma 7.4 that the power 

series expansion of Q{t) has the form 

Q{t) = qo{t - 0'+' + qi{t - 0'+'+' + ■■■ 

with qq = j^iC — qC~^)'^~^'' ■ Therefore, since go 7^ 0, the value of A when applying (4.21) 
for I{R) is X = k + I + 1. Furthermore, the gq of (4.21) is now, by [0,p.l23], equal to 

qo (C -?(-')'+' 



^pVM k\l\ «-<-i)fe+^+i- 

Thus, by applying (4.21) to I{R) and looking just at the first term of the asymptotic 
series, we obtain 

T(m - ^(-^) - r(fc + / + i)C-'^ {C-qC'f^' ^ n(T^-k-^-^\ 

' ^ y~'~ \t~vfPl + 0(i?— ), (7.13) 
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where at the last equahty we used = (. 

The asymptotics of the second integral in (7.11) is obtained by taking complex conju- 
gation in (7.13). Replacing these asymptotics into (7.11) we obtain (7.1). 

We now extend the proof to the case when u is not necessarily positive. The same 
arguments we employed above imply for tt > the slightly more general asymptotics 

Vl {VlP{-Ru + q'{x + y) + c, -Rv + q{x + y) + d)U=,, } |^^^^ = 



i fk + l 



+ o(i?-'=-'-'), 



1 



(7.14) 



for any G Q with 3|g' — g, and any integer multiples u and f of 3 at least one of which 
is non-zero. 

Suppose now v > 0. By its definition (1.3) it is clear that P{a,P) = P{P,a), so we 
can write 

P{—Ru + X + y + c, —Rv + q{x + y) + d) — P{—Rv + q{x + y) + d, —Ru -\- x -\- y -\- c) . 
Since v > 0, (7.14) applies to the right hand side above. Its asymptotics is thus 



i fk + l 



i fk + l 



c 



-l\k+l 



c 



-d+c 



k+i 



-c 



d—c 



+ o(i?-'=-'-2) 



1 



qO'+' ^ ^-d+c iC-qC-')'+' 



+ C»(i?-'=-'-2), 



1 



Jlk+l+1 



which agrees with (7.1). 

If neither u > nor v > hold, then since by hypothesis (u, v) ^ (0, 0), we must have 
—u — V > 0. By [14], P also satisfies the symmetry P{a,P) = P{—a — /3 — 1, a). Write 
therefore 

P{—Ru + X + y + c, —Rv + q{x + y) + d) 

= P{-R{-u - v) + {-1 - q){x + y) - c - d - 1, -Ru + x + y + c). 
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Relation (7.14) is applicable to the right had side above: —u — v and —v are multiples 
of 3 (since u and v are assumed to be so), (— 1 — g) — 1 = — 3 + 1 — g is divisible by 3 
(since 3|1 — q), and —u — v > 0. We obtain by (7.14) that its asymptotics is 



i fk + l 
I 



c 



i fk + l 
2^ V I 



-2c-d-l [(-l-g)C-C-T+' 

[(-«- v)C-<-i]'^+'+i 
2c+d+i [(-i-g)C-^-C]^+^ 



-C 



(1-?C)"+^ 



(« - ^;C)'=+^+l 



-c 



(w - ^;(-l)fe+«+l 



1 



+ 0(i? 



-fe-i-2\ 



-A;-Z-2\ 



(we used C ^'^ = cind — C — C ^ = 1). This is precisely the right hand side of (7.1), and 
the proof is complete. □ 



8. The evaluation of det(M") 



We evaluate the determinant of the matrix M" given by (5.11)-(5.15) by the method 
of factor exhaustion. It turns out that this evaluation holds for any C, not just for the 
specific C = e^'^*/"^ from the statement of Proposition 5.3. Therefore we will treat ^ as an 
indeterminate in Sections 8-12 and in the Appendix. 

The proof of the evaluation of det(M") involves several technical steps, presented in 
Sections 9-12 and in the Appendix. In the current section we show how to put together 
those steps to deduce the evaluation. 



^tj, and assume S >T. Then the determi- 



Theorem 8.1. Let S = YlT=i s^,T = J2 
nant of the matrix M" defined in (5.11)-(5.15), with Q, q, Xi, Hi, Wj, zj, i — 1, . . . ,m 
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J = 1, . . . ,n being indeterminates, is given by 



det(M") = (C^ - C'r'iiQ -OiQ- r^)]^-^ 

l<i<j<m 

l<i<j<n 
m n 

i=ii=i (8.1) 

Proof. Regard both sides of (8.1) as polynomials in q with coefficients in the field of 
rational functions Q{xi, . . . ,Xm,yi, ■ ■ ■ ,ym, zi, . . . , Zn,wi, . . . ,Wn, ()■ 

It is not hard to show that the degree in q of the left hand side of (8.1) is at most 
X^I^i (2*) + Sj=i (^2)- Iiideed, multiplying the first 2T columns of M" from left to right 

by 



q^ q\ q-\ q-\ . . . , q'^-\ q'^-\ q\ q\ q-\ q-\ q'^-\ q'^-\ 

q',q',q-\q-\...,q'--\q'--l^ 



(8.2) 



respectively, the degrees in q of the entries become constant along the rows (except for 
some entries generated by (5.15) when I < k). Namely, from top to bottom, these 
common degrees along the rows are 

0,0, 1, 1, . . .,Si - l,Si - 1,0,0, 1, 1, . . .,S2 - 1,S2 - 1, . . .,0,0, 1,1, . . .,s^ - l,s^ - 1, 

respectively. Thus, the degree in q of any element in the expansion of this scaled-column 
modification of M" is either or 

Sl — 1 s^ — 1 m / 

S, 



i=0 i=0 i=l 



Taking into account the factors (8.2) we multiplied the columns of M" by, it follows that 
the degree in q of the left hand side of (8.1) is at most Yl^i (2') + Sj=i (2)- 
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On the other hand, by Corollary 9.2, the left hand side of (8.1) is divisible by 

which is, up to a scalar multiple, the polynomial on the right hand side of (8.1). It 
follows that the two polynomials on the left and right in (8.1) must be the same up to 
multiplication by a factor independent of q. In particular, to prove (8.1) it suffices to 
establish it for g = 0, i.e., to show that 

det(M"|,=o) 

= ie-C'r' n [{{xi-x,)-ayi-yM{xi-x,)-c\yi-yM''' 

l<i<j<Tn 
l<i<j<n 



X 



1=1 j=l 

Taking h = in Proposition 10.1 it follows that 



(8.3) 



det(M"|,=o) = c H [{{x, - xj) - ({yi - yjMxi - xj) - c\yi - yjW' 

l<i<j<m 

l<i<j <n 
m n 

i=ij=i (8.4) 

with c G Z[C,C~^]- 

To determine c, specialize in (8.4) yi = 0, wj = 0, i — 1, . . . , m, j = 1, . . . , n. Then, 
given that we also have q = 0, formulas (5.14)-(5.15) become equivalent to 



11 



12 
21 



u 



(8.5) 
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and 



11 

12 



u 



(8.6) 



.l-k ) ^-3 _ ^3 

Denote by Mq the matrix given by (5.11)-(5.13) with definitions (8.5)-(8.6). Equal- 
ity (8.4) becomes under our speciahzation 



det(M^') = c- 



nr=in;=i(^.-^.)^ 



(8.7) 



Permute the rows of Mq by hsting them in the order 1, 3, . . . , 2S' — 1, 2, 4, . . . , 25", and 
then apply this same permutation to the columns of the resulting matrix. Denote by Mi 
the matrix obtained permuting this way the rows and columns of Mq. 

Note that when using definitions (8.5) and (8.6) in (5.12) and (5.13), in all 2 x 2 blocks 
the entries in positions (2, 1) and (2, 2) are negatives of each other. The consequence of 
this for the matrix Mi is that it has the form 



Ml 



C D 
-C C 



where C and D are S x S matrices. We obtain 

det(M^') = det(Mi) = det 



C + D D 
C 



Furthermore, (8.5) and (8.6) imply that we have 
C = iC' - ON' 

c+D^ ic' - c + r' - c^)N' = (C"' - c)(c-' + cfN', 



(8.8) 

(8.9) 
(8.10) 



where A^' is the matrix obtained from the matrix N in Theorem A.l by replacing Zj by 
-Zj, j = l,...,n. By (8.8)-(8.10) we obtain 



det(M^') = (C-^ - C)'^(C-' + C)'^(detAr')^ = (C^ - C^f^idetN'f. 
By Theorem A.l the above equality implies 



(8.11) 



det(M^') = (C^ - C"') 



2\2S 



Comparing this to (8.7) we obtain that c = (C — C ) . Replacing this in (8.4) implies 
(8.3). This completes the proof. □ 
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9. Divisibility of det(M") by the powers ofq — ^ and q — C ^ 

It turns out that rather than proving directly the divisibihty of det(M") by the factors 
involving q on the right hand side of (8.1), it is more convenient to deduce this as a special 
case of a more general result. 

More precisely, deform definitions (5.14)-(5.15) by introducing a new parameter h as 
follows: 

11 



12 

^121 



(k, I, u, v) = 



k + l 



and 



11 

12 
21 



{k,u,v) 



(9.1) 



\l-k 



- < e )a-c'f{q-Q)k,h{^-vc') 



l\l-k 



(9.2) 



where 



{a)n,h ■= a{a + h){a + 2h) ■ ■ ■ {a + {n - l)h). 
Clearly, formulas (5.14)-(5.15) correspond to the special case /i = in the above expres- 



sions. 



Proposition 9.1. Regarded as a polynomial in q with coefficients in 

cm, 

the determinant of the matrix (5.11)-(5.13) with definitions (9.1)-(9.2) admits q = 
Xh as a root of multiplicity 



^max(si - (A + 1),0) + ^max(tj - (A + 1),0), 



for all A > 0. 
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Corollary 9.2. Regard the determinant of the matrix M" given by (5.11)-(5.15) as 
a polynomial in q with coefficients in Q(a;i, . . . , Xm, yi, ■ ■ ■ , ym, zi, . . . , Zn,wi, . . . , Wn, C)- 
Then det(M") is divisible by 



Proof. Take h — in Proposition 9.1 and use that 

(m n 
max(si - (A + 1), 0) + ^ max(tj- - (A + 1), 0) 



A>0 I i=l 



J^ii^i - 1) + (s^ - 2) + ■ ■ ■ + 1] + Y^iit, _ 1) + (t^. - 2) + ■ ■ ■ + 1] 



i=l 

m 



i=l 



i=i 



□ 

Proof of Proposition 9.1. Recall the terminology we defined after the proof of Lemma 
5.2 to describe the structure of the matrix (5.11)-(5.13): the blocks outlined by large 
rectangles form block-rows and block-columns, while the 2x2 blocks form bi-rows and 
bi-columns. Refer to portion (5.12) as the ^-part and to portion (5.13) as the 5-part of 
the matrix. In block- row z, bi-row k-\- \ consists of the juxtaposition of the n + 1 matrices 



A'^^ {k,0,zj —x^ ,Wj —yi) A^^ (k ,0 , z j — X I ,zuj —yi) 



A^^ {k,tj —l,Zj —Xj^ ,vjj —yj,) A^^ (k,tj —l,Zj —x^ ,Wj —y^) 



for = 1, . . . , n, and 



B Q _'^{h,Xi^yi) B Q _rjp _^{h^Xi,yi) 



We consider first the case of the roots of type q = C ~ ^h. Perform the elementary row 
operation of adding the top half of this bi-row, multiplied by — C^, to its bottom half: 



(9.3) 



{Rp^ stands for row p of the block-row i in matrix (5.11)-(5.13)). By (9.1), the modified 
entries in bi-column I + 1 of block-column j are 



A'^^{k,l,Zj -Xi,Wj -yi)-(^A^^{k,l,Zj -Xi,Wj - yi) 

'k + l 



k 



'(.-2 _ .2 J-C-^)"+'(g-C)fc,/.(g-C)z./. 

^ ' [{zj - X,) - ({w, - y,)]'^+i+^ 



(9.4) 



and 



A^'^{k,l,Zj -Xi,Wj -yi)-C'^A^'^{k,l,Zj -Xi,Wj - yi) 



k + l 
k 



■^(i_^4)(-C-')'+^(g-C).,/.(9-C)z,/. 



(9.5) 



[{Zj - X,) - ({W^ - y^)]'^+^+' 

while by (9.2) the modified entries in bi-column Z -|- 1 of the i?-part are 

Bf\k, x., y,) - C'Bl\k, X,, y,) = - Q [({C' - OhO^Q - Ok,h{^ - K"')'"'] 
and 

Bl\k, X,, y,) - eBl\k, X,, y,) = - Q [C(l - Oi-O'iq - C)k,h{^ - K"')'"'] • 
Note that 

q-iC- Xh) I {q - Ok,h = {q-OiQ-iC-h))---{q-{C-{k- l))h) (9.6) 

provided A < k — 1. The above expressions show therefore that whenever k > A -|- 1, 
operation (9.3) produces a row along which all entries are divisible hy q — {( — Xh). 

Let D be the matrix obtained from matrix (5.11)-(5.13), (9.1)-(9.2) by applying row 
operations (9.3) for all 1 < z < m and A; > A -|- 1. Let E be the matrix obtained from D 
by performing the column operations 

C^ii C^i - C'C^i2^ 1 < J < n, / > A + 1 (9.7) 

{Cp^ stands for column p of the block-column j in matrix D). 

Then to finish the proof, it suffices to show that all entries of E along its columns 
produced by (9.7) are divisible by g — (C — Xh), and that those of them situated in bottom 
halves of bi-rows of index > A -|- 1 are in fact divisible hy {q — — Xh))'^ . 
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To verify this, consider an arbitrary entry a of a column of type (9.7) of E. 
Case 1. The entry a is in the top half of bi-row A; + 1 of block-row i. 
Then the entries of D on which (9.7) acts are unmodified by the row operations (9.3). 
Formulas (9.1)-(9.2) imply that 



a = — 



k + l 



c(i-n 



k J^'' ' ' [iZj-X,)-aWJ-y^)]''+'+^' 

Since / > A + 1, (9.6) implies that this is divisible by g — (C — \h). 

Case 2. The entry a is in the bottom half of bi-row A; + 1 of block-row i. 

Then the entries of D on which (9.7) acts are given by (9.4) and (9.5). It follows that 



a = — 



k + l 



ci-ic'-n+cii-c)] 



k J^' - ^ ' - - ^^[(^._a;.)_^(^._^.)]fe+m- 

By (9.6), this is divisible hy q — {( — Xh), since Z > A + 1. If in addition A; > A + 1, (9.6) 
implies that a is divisible by {q — {C ~ A/i))^. 

This verifies our statement about the entries of E along the columns produced by 
(9.7), and completes the proof of the statement of the Proposition for the roots of type 
q = C-Xh. 

To obtain the corresponding statement about the roots of type q = C ^ — ^h, replace 
C by C~"^ in the statement we just proved, and use the fact that this has the efi^ect of 
replacing the quantities (9.1)-(9.2) by their negatives. □ 

10. The case q = of Theorem 8.1, up to a constant multiple 

Proposition 9.1 reduces the proof of Theorem 8.1 to the case q = (see Section 8). 
It turns out that in order to handle this case it is convenient to deform the q — 
specialization of definitions (5.14)-(5.15) by introducing a new parameter h as follows: 



11 

12 
21 



(k, I, U, v) = 



k + l 

k 




1 



{u - V()k+l+l,h 




{U-VC ^)k+l+l,h 



and 



B, 
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{k,u,v) = 



(10.1) 




l — k,h 



(10.2) 



where as in the previous section {a)n,h '■= ci{a + h){a + 2h) ■ ■ ■ {a + {n — l)h). Clearly, 
setting h = above yields the specialization g = of the definitions (5.14)-(5.15). 

Proposition 10.1. Let Mq be the 2S x 2S matrix given by (5.11)-(5.13) and (10.1)- 
(10.2). Then we have 

det(Mo) = 

8i-lSj-l 

^ n Yl Yl[i^i-^3)-Ciyi-yj)-ik-^)h][{xi-Xj)-C-^{yi-yj)-{k-l)h] 

'i-<i<j<m k=0 1=0 
ti-ltj-1 

^ n n n t"^^* ~ ^-j) ~ ~ "^j) - - o^][(^i - - c^i^i - -ik- i)h] 

l<i<j<n fe=0 1=0 
m n Si — ltj — 1 

^nn n n {[(^i-^.)-c(yi-«^,)-(/^+o/^]-' 

i=lj=l k=0 1=0 

x[{xi - Zj) - C\yi - Wj) -{k + l)h]-^] , (10.3) 

where c G C"-*-] . 
Proof. Let 

m n Si+ij— 2 

c^ = nn n [{xi-Zj)-C,{yi-Wj)-ah][{xi-Zj)-C^{yi-Wj)-ah]. (10.4) 

i=lj=l a=0 

It follows from the definition of Mq that when multiplying it by d, all its entries become 
polynomials in the variables h,xi, . . . , Xm, Hi, ■ ■ ■ , Vrm -^i, ■ • • , wi, . . . ,Wn with coefii- 
cients in Z[C, C""*"]- 

For a e Z and integers s,t>l define 

Ni{a,s,t) = \{{k,l) : k - I = a,0 < k < s - 1,0 < I < t - 1}\ (10.5) 

and 

N2{a, s, t) = \{{k, I) : k + I ^ a,0 < k < s - 1,0 < I < t - (10.6) 
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Then, multiplying all rows of Mq by cZ, (10.3) becomes equivalent to 
det(cZMo) 

I 

^<i<3<rn a=—{sj — l) 
l<i<j<n a=-{tj-l) 

xii^i - zj) - C\wi - Wj) - a/i]^^(«'**'*^)} 

m n Si+tj—2 

xiin n {[ix^-ZJ)-ay^-w,)-ahr-''^(^'^-'^^ 

i=lj=l a=0 

(10.7) 

Regard both sides of (10.7) as polynomials in the variables 

hi Xi, . . . , Xmj yii ■ • • 1 ymi Zi, • • • , Zfii ■ ■ ■ i '^nj 

with coefficients in Z[C,C~^]. 

By Corollary 11.3 and Proposition 12.1, det((iMo) is divisible by the full product 
expression on the right hand side of (10.7). It follows from definitions (10.5) and (10.6) 
that 

s-l s+t-2 

Ni{a,s,t)= J2 N2{a,s,t)^st. 

a=-{t-l) a=0 

By (10.4), the degree of d is 

m n 
i=l 3 = 1 

Therefore, jointly in its variables, the polynomial on the right hand side of (10.7) has 
degree 

m n 

2Sdeg{d) + 2 ^ SiSj + 2 ^ Utj - 2 ^ Sitj. (10.8) 
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To complete the proof, it suffices to show that the degree of dct(dMo) is less or equal 
than the value given by this expression. In turn, this is clearly equivalent to the degree^ 
of det(Mo) being at most the value given by the three sum-terms in (10.8). However, 
this follows by the argument we used in the proof of Proposition 5.3 to factor out the 
powers of R from the determinant of the matrix with blocks (5.16)-(5.17): Indeed, the 
degree in R of each entry of that matrix is equal to the joint degree in its variables of the 
corresponding entry of Mq. □ 

Remark 10.2. Setting /i = in (10.3) and using (8.3) (which is a special case of Theorem 
8.1) we obtain the explicit value c = (C~^ ~ C^)^'^ fo^^ the multiplicative constant in 
Proposition 10.1. This provides a more general deformation of Cauchy's determinant 
than the one presented in Theorem A.l. 



11. Divisibility of det{dMo) by the powers of {xi — Xj) — C^^iVi ~ Vj) ~ 

and {zi — Wj) — C^^{zi — Wj) — ah 

For any integers s > 1 and < i < s — 1, define the 2s- vector v| by 
v| = Ko, -C^Vio, Vil, -C^va, . . . , Vs-i, -C^Vj,s-i], 

Vij = i-iyjl Q h^, j = 0, 1, . . . , s - 1. (11.1) 

Denote by 0,- the r- vector with all coordinates equal to 0. 

Lemma 11.1. Let k and I he integers, 0<k<Si — l,0<l<Sj — l,i<j. Then, if 
+ CiUi " Vj) -\- {k — l)h, the linear combination of the rows of matrix Mq specified 
by the 2S-vector 

[02si, • ■ ■ : 02s,_i, -V^% 02si+i, • • • , 02sj_i, V;""', 02s^. + i, " " " , 02s^] (11-2) 

is a vanishing row combination. Furthermore, for any fixed a e Z, the vectors of the 
above type satisfying k — I = a are linearly independent. 

Proof. Note that in any subset of the set of vectors with k — I = a fixed there is 
a unique vector with a rightmost non-zero coordinate. Indeed, if two such vectors have 

^The degree of a rational function is the difference between the degrees of its numerator and deno- 
minator. 
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their rightmost non-zero coordinates in the same position, the values of I they correspond 
to are equal, but then k — I = a implies that the /c-values they correspond to are also 
equal, so the two vectors are identical. This implies the linear independence claimed in 
the statement of the Lemma. 

Due to the form of Mq, the statement that the specified row combination is vanishing 
amounts to four separate identities. Namely, wc need to check that the indicated row 
combination produces along the even and odd indexed columns of portion (5.12) of 
Mo, and also along the even and odd indexed columns of portion (5.13) of Mq. It is 
convenient to refer to these two portions as the yl-part and S-part of Mq. 

Along column 2r + 1 of block-column h (cf. our terminology related to the structure 
of (5.12)-(5.13)) of the ^-part of Mq, the vanishing of the indicated row combination 
amounts to 

Si-l 

- ^Vka[A^^{oi-, r, Zb-Xi,Wb- Vi) - C^^^^(q!, r, Zb - xi, Wb - Vi)] 

Sj-1 

+ XI ^«a[-4^^(a,r,2;6 - xj.Wb - yj) - C'^A'^^{a,r,Zb - Xj,Wb-yj)] = 0, 

a=0 (11.3) 



provided Xi = xj + ({yi — yj) + (k — l)h. 

By (11.1) and (10.1)-(10.2), the left hand side above equals 

Ti k\ /g + ry-C 1 

^ {k-a)\\ a J C {(zb-Xi) -C{wb-yi))a+r+i,h 

+ Yi-hr-^^f'' + 1^ . 

^ {l-ay.\ a J C {{zb - Xj) - C{wb - yj))a+r+i,h (^^^^ 

When Xi = Xj + ({yi - yj) + {k - l)h, we have {zb - xj) - ({wb - yj) = {zb - Xi) - ({wb - 
yi) -\- {k — l)h. Denoting c—{zb — Xi) — C{wb — yi), we obtain that (11-3) is equivalent to 

^ k\ f a + r\ 1 



^0 {k-ay.\ a J{c)a+r+i,h 



^0 (^-«)!V « J ic+ik-l)h)a+r+l,h 



(11.5) 
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We claim that the left hand side of (11.5) evaluates as 

k 



k\ ( a.^r 



1 



1 



a=0 {k-a)\\ a y (c)a+r+l,/i h'^^{^^k)r+\ 

Indeed, let c' = cjh. The left hand side of (11.6) can then be written as 

k 



(11.6) 



{k-Oi)\\ a ; (c')«+r+l/i'*+^+^ 

fc 



- Xl( ^)«*^ 



1 



a=0 
1 



ct! (c')r+i(c' + r + l)a 
— /c, r + 1 



c' +r + l 



;1 



Then (11.6) follows by Lemma 11.2. 

By (11.6), the right hand side of (11.5) equals 



1 



1 



This proves (11.5), and therefore (11.3). 

Next, we verify that the row combination (11.2) is vanishing along the column 2r + 2 
of block-column b. This amounts to proving a close analog of (11.3); namely, that 

Si — l 



a=0 



^Vka[A^'^{a,r, Zb-Xi,Wb- yi) - C^A^^(q;, r, Zb - Xj, Wb - Vi)] 

8j-l 

+ XI ^«a[^^^(a, r, Zb - Xj,Wb -Vj)- (^A^'^{a, r, Zb - Xj,Wb - Vj)] = 0, 



a=0 



(11.7) 



provided Xi — Xj + ({i/i — pj) + {k — l)h. By (11.1) and (10.1)-(10.2), the left hand side 
of (11.7) is nearly the same as the left hand side of (11.3) — the only difference is that the 
factor 1 — at the numerator gets now replaced by — C^- Therefore (11.7) follows 
by (11.3). 

To complete the proof we need to check that row combination (11.2) vanishes along 
the columns of the S-part of Mq. 
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Along column 2r + 1 of the S-part of Mq, the vanishing of row combination (11.2) 
amounts to showing that for Xi = Xj + ({yi — yj) + {k — l)h one has 



S;-l 



a=0 



^ L'A;a[-B"(Q!,r, Xi,yi) - C'^B^'^{a,r,Xi,yi)] 

Sj-l 

+ ^ via[B^^{a,r,Xj,yj) - C^B^^{a,r,Xj,yj)] = 0. 



a=0 



By (11.1) and formulas (10.1)-(10.2), the left hand side of (11.8) becomes 



Letting c = (xi — Cyi)/h, the first sum above can be written as 



a=0 



{k — a)\ \a 



a=0 



(c)./i"2Fl 



(-r)a (c)r 

a! (-c-r+ 1)q 

— /c, — r 
-c - r + 1 ' 



(11.8) 



(11.9) 



Replacing r by — r — 1 and c by — c + 1, Lemma 11.2 implies, after simplifications, that 



—k, —r 



-c - r + 1 ' 



;1 



(c - k), 

(c)r 



Therefore, the first sum in (11.9) equals h'^{c — k)r- But the second sum in (11.9) has the 

same form, except {xj — C,yj)/h = {xi — Cyi)/h — {k — l) = c— (A; — /), since we are assuming 
Xi = Xj +C{yi — yj) + ik — l)h. It follows that the second sum equals h'^ {{c — {k — I)) — 
and thus cancels the first sum. This proves (11.8). 
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The final identity to check is that row combination (11.2) is vanishing along column 
2r + 2 of the S-part of Mq. This is equivalent to the following counterpart of (11.8): 



Si-l 



a=0 



Sj-l 

+ XI via[B^'^{a,r,Xj,yj) - C'^B^^{a,r,Xj,yj)] = 0. 



(11.10) 



a=0 



Replacing (11.1) and (10.1)-(10.2) in this shows that, just as it was the case for (11.3) 
and (11.7), the expressions on the left hand sides of (11.10) and (11.8) differ only in that 
the factor 1 — in the former gets replaced by — the latter. Therefore (11.10) 
follows from (11.8), and the proof is complete. □ 

Lemma 11.2. Let < A; e Z, r e Z. Then if c is an indeterminate one has 



2^1 



—k, r + 1 
c + r + 1 



(c) 



r+l 



if r + 1 > 0, 



{c+k)r+l' 

(c + r + l)_,_i ' it^ + -L<U- 



Proof. Since k > 0, the hypergeometric series is terminating and the equality in the 
statement is equivalent to a polynomial identity in c. By Gauss' formula (4.9), for any 
large enough integer c we have 



-k,r+l^ 
c+r+1 ' 



r(c + r + l)r(c+A;) 
r{c)T{c + k + r + l)' 



The expression on the right hand side readily checks to simplify to the ones in the 
statement. This proves that the rational functions in c on the two sides of the equality in 
the statement agree on infinitely many values of c. They must therefore be identical. □ 

Corollary 11.3. Regarded as a polynomial with coefficients in Z[(, (~^] in the variables 
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h, xi, . . . ,Xm,yi, ■ ■ ■ ,ym, zi, . . . , Zn,wi, . . . ,Wn, det{dMo) IS divisible by 

n ii {[(^^ - - - yj) - ahf^^-'^-^^^ 



^<i<3<:m a= — (Sj — 1) 



ti-1 



x[{x, - xj) - c\yi - yj) - a/i]^^('^'*^'^^)} 



X n n {[{zi - zj) - Ciwi - wj) - ahf^^-''-'^^ 

l<i<j<n a=-{tj-l) 

x[{zi - Zj) - C\w, - Wj) - a/i]^i('^'**'*^)} . 

(11.11) 

Proof. By Lemma 11.1, Xi = Xj + ({yi — yj) + ah is d: root of multiplicity A^i (a, Sj, sj) of 
det((iMo), when the latter is regarded as a polynomial in the single variable Xi, with co- 
efficients in Z[CX~^^h,Xi,...,Xi-i,Xi+i,...,Xm,yi,---,yTn,Zi,...,Zn,Wi,...,Wn\. By 

Bezout's theorem, it follows that 

[{xi - xj) - Civ, - Vj) - a/i]^^("'^-^^) I det(rfMo). (11.12) 

By the form of the expressions (10.1)-(10.2) and (10.4), the effect of replacing C by 

in dMo is to replace each entry by its negative. Therefore (11.12) implies that we also 

have 

[{xi - Xj) - C\yi - yj) - a/i]^^^'^'^^'^^) I det(dMo). 

To obtain the divisibility by the factors on the second line of (11.11), note that Lemma 
11.1 admits a perfect analog for column combinations. The proof of this analog is identical 
to the part of the proof of Lemma 11.1 that pertains to verifying that combination (11.2) 
is vanishing along the columns of the A-part of Mq. This is due to the symmetry of 
formula (10.1) under interchanging k and I. □ 

12. Divisibility of det(ciMo) by the powers of {xi — Zj) — C^^(j/i — Wj) 

Proposition 12.1. Regarded as a polynomial with coefficients in Z[C,C~^] in the vari- 
ables h, xi, . . . , Xm, yi, ■ ■ ■ , y-mi zi, . . . , Zn, Wi, . . . , Wn, det((iMo) is divisible by 

m n Si+tj—2 

nn n K^^-^o-ai/.-^^.o-a/i]^^-^^^"'^-*^) 

i=l j=l 0=0 

X[(X, - Zj) - C\y^ - W,) - a/z]25-iV.(a,.„t,)_ 

61 



Our proof is based on finding suitable sets of independent finear combinations of the 
rows and columns of dMo that vanish when {xi — Zj) — C^^iVi — Wj) — ah = 0. These are 
presented in the following four preliminary results. 

Lemma 12.2. Each column of dMo outside of block- column j becomes zero when {xi — 
Zj) - C^^iVi - Wj) - a/i = 0. 

Proof. This holds on account of the multiplicative factor d — outside block-column j 
the factors {xi — Zj) — C^^(2/i ~ ^j) ~ '^h of d do not get canceled by factors in the 
denominators of expressions (10.1)-(10.2). □ 

Lemma 12.3. Let a e {0, 1, . . . , + — 2} be fixed, and assume {xi — Zj) — C^^(2/i ~ 
Wj) — ah = 0. Then the column combination 

[02ti, • • • ) ^2tj-i, O2I, 1, — C^^) 02tj-2l-2, 02tj+i, • • • ) 02t„, 02S-2t] 

in dMo is vanishing, for I = 0, . . . ,tj — 1. 

Proof. The vanishing of the given column combination along row 2r + 1 of block-row 
i of dMo amounts to 

d''A^\r, I, Zj - Xi,Wj - Vi) - C^^d''A^\r, I, Zj - Xi, Wj - pi) = 0, (12.1) 

where 

d'^ = n [{xi - Zj) - CiVi - Wj) - ah][{xi - Zj) - C\yi - Wj) - ah] (12.2) 
0=0 

is the part of d involving precisely indices i and j. By definition (10.1), (12.1) is equivalent 
to 

r + l 
r 



X 



^-1(1 _ ^±2^-2x ii^j - ^i) - C(Wj - y^))s,+t,-lM{i^J ~ ^t) ' C H^j - |/i))s,+t,-l,fe 

{{Zj - Xi) - CiWj - yi))r+l+l,h 

_ ^/^ _ ^±2^2^ ii^j - ^i) - CK- - y^))s,+t,-l,h{i^j - ^i) - C~\wj - yi))s^+t^-l,h 

{{Zj - Xi) - C~^{Wj - yi))r+l+l,h 

= 0. 
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(12.3) 



Suppose +2 is chosen at the exponent in (12.3). Then the first term in the square brackets 



is clearly 0. The second term is due to the factor of {zj 



C{wj — yi) + ah in the 



numerator, which clearly does not get canceled by any factor of the denominator. The 
case when —2 is chosen at the exponent in (12.3) follows similarly. 

The vanishing of the indicated column combination along row 2r + 2 of block-row i of 
dMo amounts to 



d^^A^^{r, I, Zj - Xi, Wj - Ui) - C d'-'A^^(r, I, zj - Xi, wj 
This is obtained from (12.1) by dividing by — C^^- 



Vi) 



0. 



□ 



Lemma 12.4. Let a — Si + k, k e {0, . . . ,tj~2}, be fixed. Then when {xi — Zj)—(^^{y i — 
Wj ) — ah = 0, columns 2, 4, 6, . . . , 2/c + 2 of block-column j of dMo become identically zero. 

Proof. For < I < k and < r < — 1 , the entry in column 2Z + 2 of block-column j 
of dMo that is situated in row 2r -|- 1 of block-row i is 

d'-^A^^ir, I, Zj - Xi, Wj - yi), 

where d'^^ is given by (12.2). By (10.1) and (12.2) the explicit value of this is 



r + l 
r 



-3 ii^j - ^i) - C{wj - yi))si+t^-i,h{{zj -Xi)-C ^{wj - yi))si+t^-i,h 



{{Zj - Xi) - C{Wj - yi))r+l+l,h 



{{zj - Xi) -({wj -y^))s. 



+ti 



X 



:) -C ^{u!j - yi)) 



Si+t-j —l,h 



{{Zj -Xi)-C ^{Wj - yi))r+l+l,h 



ah in the numerators are not canceled by any of the 

■ hh with 



■Wi 



The factors {xi — Zj) — C, {y 

factors in the denominators, as the latter are of the form (xi — Zj)— iUi ' 
< b < r+l < Si+k — a. Thus the expression above is for {xi—Zj) — ('^^{yi—Wj) — ah — 
0. 

The entry in column 2/ -|- 2 of block-column j of dMo which is situated in row 2r + 2 
of block- row i is seen to be by the same argument. □ 

Lemma 12.5. Let a e {0, . . . ,tj — 2} be fixed. For 0<k<tj — 2 — a, consider the vectors 



{-i)Htj-i)i 

k\{tj-2-a-k)\ 



{-hy 



and 



Vfc 



02fe,0, 



(t,-l)! 



{a — a)\{a — a -\- 1 -\- k) 



— 1^) ) "2(tj-2-a-fc) 



0<a<a 



{tj-2-a- k)\ 
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Then if {xi — Zj) — ^^^{yi — Wj) — ah = 0, the column combinations ofdMo with coefficients 

[02ti 02t,-i , Ufc, Vfc, 02t,.+i 02t„] (12.4) 

are vanishing for all k = 0,1, ... ,tj — 2 — a. 

Proof. The fact that the indicated column combination is vanishing along row 2r + 1 
in block-row i is equivalent to 

(-l)^(t,-l)! ^ ^,12, 
k\{tj -2-a-k)\ ^^{a-ay.{k+l + a-a) ^' j y*/ 

= ,,, {-hr+^^HA^\r, a + k+l,z^- x^w^ - y,), 

[tj — z — a — K). (12.5) 

provided {xi — Zj) — C^^iVi ~ — ah = 0. 

Assume {xi—Zj)—({yi—Wj)—ah = 0. Using (10.1), (10.4), the fact that {a)n,h/ {o)k,h = 
{a + kh)n-k,h and dividing through by all the factors of d not present in (12.2), the above 
equality becomes 

(-l)"-(^j - 1)! (~/i)" /r + a' 

k\{tj - 2- a- ky. ^^{a- a)l{k + l + a- a)\ a 

X [C~^{-ah +{r + a + l)h)si+tj-r-a,h{u - vC~^)si+t,+i,/i 

-C^{-ah)s^+t^+i,h{u - vC^ + {r + a + l)h)s^+t^-r-a,h] 

^ _^(!iIll)l^r_Ma+fc+i fr + a + k + 1 
{tj-2-a-k)V ' V ^ 

X [C^i-o-h + (r + a + A; + 2)h)si+tj-r-a-k-i,h{u - vC~^)si+t,+i,h 

-C^{-ah)s^+ti+l,h{.u - vC^ + (r + a + /c + 2)h)si+tj-r-a-k-i,h] , 

where u = zj — Xi and v = wj — yi. Using that {—ah)si-\-t -\-i,h = 0, the above simplifies 
to 



^ (a-a)!(/c + l + a-a) V « / ' 



= (-l)U! f"" + ^ + ^ + ^Vr + /c + 2),,+,^._,_„_fe_i. 

^ ^ (12.6) 
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Write for notational simplicity s = Si + tj. Then after some manipulation, the left hand 
side of (12.6) becomes 

{s - ay. 'A (-a)a(-Q -k - l)a{r + 1)q (r + 1 - a)a 

a\T{r + l){-a-k-l) a\{-a-k)a (r + 1 - a)^ ' 

an expression that is well-defined for any 0<a, /ceZ, a<seZ and r G C. Furthermore, 
when r ^ Z, the above expression can be written as 

{s-a)\{r + l-a)a -A (-a)a(-a - k - l)^{r + 1)^ 
aW(r + l)(-a - k - 1) ^ a\i-a - k)Jr + 1 - a)^ 



(g - a)!(r + 1 - a)a 
' a^T{r + l){-a-k-l) 



3-^2 



—a, —a — k — l,r + 1 
r + 1 — a, —a — k 



(12.7) 



The 3F2 above can be evaluated by choosing n = a, x = —a — k — 1, y = r + l, z = r+l — a 
in the Pfaff-Saalschutz summation formula (found e.g. in [23, (2.3.1.3), Appendix (III. 2)]) 



3^2 



-n,x,y 
z,l + X + y — z 



n 



{z - x)n{z - y), 



{z)n{z - X - y)n 

Thus, by (12.7) we obtain that for r ^ Z the left hand side of (12.6) equals 

{s-a)\{r + l-a)a {r + k + 2)a{-a)a 



a!r(r + l)(-a - /c - 1) (r + 1 - a)a{k + l)a 

However, this is readily seen to agree with the right hand side of (12.6) (which, by writing 
ix + 'n)m-r = r(w + n)/r(r + n), is defined for any r G C \ Z, m, n G Z). This proves 
equality (12.6) for r G C\Z. Since the two sides of (12.6) are analytic functions of r that 
coincide on C \ Z, they must be identical. This proves (12.6) and therefore (12.5). 

The vanishing of column combination (12.4) along row 2r+2 in block-row i is equivalent 
to the equality obtained from (12.5) by replacing the A^^'s by A^^'s. This is readily seen 
to be equivalent to (12.6) as well (the only difference from the case of A^^'s is that 
the factors and in the formula displayed above (12.6) get replaced by and 
C, respectively). This completes the proof of the statement of the lemma for the case 
{xi - Zj) - CiUi - Wj) - ah = 0. 

The case when {xi — zj) — C~^{yi ~ Wj) — ah = follows by replacing ^ by in 
the case we proved above, and using the fact that the right hand side of (10.1) simply 
changes sign under this substitution. □ 
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Proof of Proposition 12.1. We treat first the 

For each a G {0, . . . , si + tj — 2}, Lemmas 12.2 and 12.3 provide together a set of 
25' — 2tj + tj = 2S — tj independent column combinations that are vanishing when 
{xi — Zj) — C^^iUi — Wj) — ah = 0. To complete the proof it suffices to enlarge this set by 
tj — N2{a, Si, tj) more vanishing column combinations, so that the enlarged set of column 
combinations is linearly independent. 

If < a < — 2, Lemma 12.5 provides tj — 1 — a such additional column combinations. 
However, it follows by definition (10.6) that N2{a, Si,tj) = a + 1 for < a < tj — 2. 
Therefore these can be taken as the additional tj — N2{a,Si,tj) column combinations 
that we need. 

For tj — 1 < a < Si — 1, (10.6) implies that N2{a, Si,tj) = tj, so no vanishing column 
combinations are needed in addition to the ones supplied by Lemmas 12.2 and 12.3. 

If on the other hand Si < a < Si + tj — 2, then N2{a, Sj, tj) — Si + tj — a — 1, and the 
needed additional vanishing column combinations are provided by Lemma 12.4. 

Similar arguments, based on the counterparts of Lemmas 12.2-12.5 presented below, 
imply the same conclusion if Si <tj. □ 

Lemma 12.2'. Each row of (IMq outside of block-row i becomes zero when {xi — Zj) — 
C^^iVi - Wj) -ah = 0. 

Proof. This follows by the same argument as Lemma 12.2 □ 

Lemma 12.3'. Let a e {0, 1, . . . , Si + tj — 2} be fixed, and assume {xi — Zj) — C^^iVi ~ 
Wj) — ah = 0. Then the row combination 

[02si, • • • ) 02si_i, 02r, 1, -(^'^ ,02si-2r-2,02si+i, ■ ■ ■,02sm] 

in dMo is vanishing, for r = 0, . . . , Sj — 1. 

Proof. The vanishing of the given column combination along column 21 + 1 of block- row 
i of dMo amounts to 

d''^A^^{r, I, Zj - Xi, Wj - yi) - (^"^d^^ A^^ {r, I, Zj - Xi, Wj - yi) = 0. 

with (i*-' given by (12.2). This follows by the argument used in the proof of Lemma 12.3. 
The case of column 21 + 2 follows by dividing the above equality by C^^. 

In addition, we need to verify the vanishing of the indicated row combination for the 
S-part of dMo- But this holds trivially since the entire S-part of cZMq is zero when 

{xi - Zj) - C^^{yi - Wj) - ah = 0. □ 

Lemma 12.4'. Let a — tj+l, I E {0, . . . , — 2}, be fixed. Then when {xi — Zj) —C^^iVi ~ 
Wj) — ah = 0, rows 2, 4, 6, .... 2/ + 2 of block-row i of dMo become identically zero. 

Proof. This follows by an argument perfectly analogous to the one used in the proof 
of Lemma 12.4. □ 
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Lemma 12.5'. Let a e {0, . . . , Sj — 2} be fixed. For 0<k<Si — 2 — a, consider the 
vectors 



k\{si - 2 - a - ky. 



0, 



i-hy 



{a — ay.{a — a + 1 + k) 



0<a<a 



and 



02fc,0, 



{si-2-a-ky. 



Then if (xi — Zj) — C (Ui — wj) — ah = 0, the row combinations of dMo with coefficients 

[02si, . • • ) 02si_i,Ufe, Vfc, 02si+i, • • • , 02sm] 



are vanishing for all k = 0,1, Si — 2 — a. 

Proof. The vanishing of the indicated row combination on the ^-part of dMo follows 
by the calculations in the proof of Lemma 12.5. Along the B part of dMg this holds 
trivially, as explained in the proof of Lemma 12.3'. □ 



13. The proofs of Theorem 2.1 and Proposition 2.2 

We are now ready to present the proofs of the two results from which we deduced 
Theorem 1.2 in Section 2. 

Proof of Theorem 2.1. By reflecting the lattice and the holes about a vertical lattice 
line if necessary, we may assume S = YlT=i ^ X^/Li U = T. 

The statement follows then directly by Proposition 5.1, Proposition 5.3 and Theorem 
8.1, using that for C = e^^*/^ one has 

r'-c' = ^\/3 (13.1) 

and 

{u-Cv){u-C^v) = u'^+uv + v^. (13.2) 

□ 

Proof of Proposition 2.2. By Proposition 5.1 we have 



Cj{E{ai,qai), . . . ,E{as,qas)) = |det(Mi)|, 
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(13.3) 



where Mi is the matrix 



Uo{ai,qai) i7o(oi -l,goi + l) 
;7o(ai + l,gai-l) J7o(ai,qai) 



Us-i(ai,qai) (oi -l,qai + l) 

C/s_i(ai + l,gai-l) Us-i(ai,qai) 



Uoias,qas) Uo{as — l,qas + l) • ■ • Us-i{as,qas) Us-i{as — l,qas+l) 

.Uoias + l,qas — l) Uo{as,qas) ■■• Us-i{as + l,qas — l) f^-i (a^ .qa^) 

Apply the determinant-preserving operator (5.4) to Mi twice, first along the rows of 
odd index and then along the rows of even index. By the arguments in the proof of 
Proposition 5.3 we obtain that 



det(Mi)= Yl {ak-aifdet{M2), 

l<k<l<s 



(13.4) 



where M2 is the matrix 



'D"Uo(aj^,qai) 
■D°Uo(ai+l,qai-l) 



■D°Uo(ai-l,qai + l) 
■D°Uo(ai,qai) 



'D^U^_l(ai,qa-i) 
■D°Us_l(ai+l,qai-l) 



C°C/s_i(ai-l,qai + l) 
T>°Us-l(ai,qai) 



We have seen in the proof of Proposition 6.1 that 

Ui{a,b) = - ^[C-'-\a - bCY - C^'+'ia - bC'Y] 

+ monomials in a and b of joint degree < /. 

This implies that, viewed as a block-matrix consisting of 2 x 2 blocks, M2 is block-upper- 
triangular. Furthermore, using Lemma 6.4 and the assumption 3\1 — q, the {k + l)-st 
diagonal 2x2 block of M2 is seen to be 

■^-1(1 _ qQk _ _ q^-i^k ^-3(1 _ q^^k _ ^3(1 _ q^-y 

.c(i - ^0' - r^i - c-'(i - qO" - C(l - qc-')' . 

Since the determinant of this block equals 

1 



^ 

'2^ 



47r2 



ic'-cnq-cnq-cr, 



equality (2.2) follows by (13.3) and (13.4), using (13.1)-(13.2). 
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□ 



14. The case of arbitrary slopes 

Let us now consider the pure linear multiholes of not necessarily equal slopes , • • • , 

and VF^^, . . . , W^". Then their correlation is still given by Lemma 5.1, but with the fol- 
lowing change in the definition of the block matrices Cj and D of matrix M: Modify the 
blocks of (5.2) so that the top left entries of the 2x2 matrices change from 

P(-R^.-■^j-^i) + °■i\-bj\-\,-R(wj-y^)+q(a^'^-bj-^)-\) ••• P {- R{z ^ - x ^) ^cn^-b ^ ^ ^ - \ ,- R^-w j - y i) + q{a^x-^ j ,t ^) 

P(-R(zj -Xi) + ai2-bji-l,-R(iVj - yi)+q(ai2-bjl)-l) ■•• P(-R(zj -Xi) + ai2-bj^tj -Hi) + 9(ai2-bj,t^. )-l) 



P(,-R(zj -a:i) + ai g^ -b^i -l,-Jl{-u;j -3/j)+q(oj_s^ -bji)-l) 

to 

P{-R{zj -x^) + a^-J^-bJ-J^-l,-R(nlj - y ^)+q^a^'^^- q'_.b -1) 
P(-R(zj -x^)+ai2-bji-l,-R(iVj-yi)+q^ai2-'ljbji-l) 



P(-R(zj -Xj) + aj - -a{) + 9(«i,Sj -bj,tj)-'^) 

P(-R{zj -aJjj + aji-bj' t^. -1,-R(nij - S/j ) + a j j - q^. b^' j^. -1) 
P(-R(zj -Xi) + ai2-bj^tj - ^'-^i'>"j-yi) + 1i°-i2-l'jbj^t^ -1) 



P{-R{zj -1:^) + a^g^ -6ji-l,-Jl(«)j-3/j)+qjaj -q^. 6^1-1) ••• P{-R{zj -x^)+a^,|^ ~''->~^('^j ~yi'>+1i°'i,Si -1jbj,tj -1) 

and the blocks of (5.3) so that the top left entries of the 2x2 matrices change from 

Uo(R^i+aii,Ryj+qaii) ■■■ Ug _i(Rxi+aii , Ryj +qaii) 

Uo(Rxi+ai2,Ryj+Qai2) U^S-T- 1 (^^i +«i2 . +9ai2) 

to 

UoiR^i+a-U'Ryj+li^-il) ••• Us-T-li^^i+^-il'^yj+Qi^-il) 

UQ(Rxi+ai2,Ryj+<!ia,i2) ••' C^S-T-l (fl^i+'^i2 .^fj +'Jiaj2) 



UQ(,Rx^+a^ ,,^ ,Ryj+q^ai g^) 



Us-T-l(l^^i+°'i,Si'l^yj+''i°'i,Si) 
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It is apparent from the above that in this more general version of the matrix M wc have 
the same 2m + 2n opportunities to apply operation (5.4) (or its analog for columns) as 
for the original matrix M of Section 5, obtained from the above when qi — ■ ■ ■ — q'n — Q- 
Indeed, the arguments of Section 5 apply in this more general case of arbitrary slopes, 
and give that the determinant of the above modification of the matrix M is equal to the 
determinant of the deformation of the matrix M' of (5.6)-(5.9) in which the quantities 
qaik are replaced by qittikS and the quantities qbji by q'^hjis. The asymptotics of the 
thus modified entries in (5.9) is given by Proposition 6.1. For the modified entries of 
(5.6) we need the following slight extension of Proposition 7.1, whose proof follows by 
the arguments of Section 7. 

Proposition 7.1'. Let q and q' he rational numbers so that 3|1 — q and 3|1 — q' . Then 
for any fixed c, d & I^, u,v & 3Z, {u, v) ^ (0, 0), and < A;, Z e Z, we have 

V[ {vlP{-Ru + x + y + c,-Rv + qx + q'y^d%=a,]\y^^,^ = 



~2^\ I 



:-d-i {l-qCr{l-q'Cy 
{u - ^;C)fc+«+l 



1 



The arguments that proved Proposition 5.3 yield then that 

uj{El\{Rxu Ryi), . . . , El2{Rxm,Rym).W^l{Rzi, Rw^), W^-jRzr,, Rw^)) 

= {h) i[ n K-«-)^n n c 



where 



i=l l<j<fc<Si 1=1 \<j<k<ti 

X |det(-M")| i?^'f^l<i<J<r« *<*J+5Zl<i<j<n Zli^l Z]j = l 

M" = \C" D"], 



(14.1) 
(14.2) 



C" 



(14.3) 
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r" 



Al-(0,0,Zj-a:i,wj-Vi) 



A^^(0,0,Zj-Xi,Wj-yi) 



{0:tj-l:Zj-Xi,wj -Vi) Alj (0,tj-l,Zj -Xi ,^ ^ - 



('^i-^'0,ZA-Xi,WA-yi) A (s j - 1 ,0,z,- -a: j -j/j ) 



^i,- {si-'^,0,Zj-xi,WA-yi) A.. (si-l,0,Zj-xi,iVj-Vi) 



Alj(.Si-l,tj-l,Zj-Xi,-Wj-yi) A\'^(si-l,tj-l,Zj-Xi,Wj 



Alj{si-\,tj-l,Zj-Xi,Wj-yi) A\l(ai-l,tj-l,Zj-Xi,Wj 



D" 



(14.5) 



and 



0,i 



0,i 
,21 



As^-l,x^,y^) 





11 


< 


12 


A. 


21 


13 



,12 . 
,22 , 



.(Sj-l,xj,i;j) 







[1 









^l^.-r-i.iCo.^i.fi) 



_1 j(0,Xj,!/j) 



Vi) 



-1 i(^i-l>^i>2'i) - 



(14.6) 



_ji-Q,c)'(i-?K)' 



(w - 

c 



(W - -i;(-l)fe+«+l 



(14.7) 
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11 



12 
21 



-l\l-k 



(14i 



Theorem 3.1 — the asymptotics of the correlation in the case of arbitrary slopes — follows 
by (14.1) and the following extension of Theorem 8.1. 

Theorem 14.1. Let S = J^iLi ^i, T = YJj=i tj> (^i^d assume S >T. Then the determi- 
nant of the matrix M" defined in (14.2)-(14.8), with q, Xi, yi, wj, zj, i = 1, . . . ,m, 
j = 1, . . . ,n being indeterminates, is given by 

m n 

det(M") = (c^ - c-^)^^ ntfe - ofe - c-^)]^^^) n[(^^- - ow, - c')p^ 

i=i j=i 

l<i<j<m 
l<i<j<n 



n Hid^^ - ^j) - C{yi - Wj)){{xi - Zj) - c \yi - Wj))]- 



-Siti 



i=lj=l 



(14.9) 



Proof. Theorem 8.1 was proved in Sections 8-12. The only part of the proof affected 
by the generalization from the matrix M" to M" is Section 9, as we now need to show 



that detM" is divisible by {qi — C )^^'^ and {q'j — C )^^'^- As in Section 9, deform the 
definitions (14.7) and (14.8), by introducing a new parameter h, to 



11 

12 



. V 21 J , , fk + 1 

{k,l,u,v) = ( ^ 



















i] 
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{u - vC)''+i+^ 

{U - -f;(-l)fc+«+l 



(14.10) 



and 




{k,u,v) = 



l-k 



- { e )a-C'f{Qr-OkAu-vC') 



l-k 



(14.11) 



where {a)n,h '■= o(a + h){a + 2h) ■ ■ ■ (a + (n — l)h). Running the arguments of the proof 
of Proposition 9.1 with the sole change that the row operations (9.3) and (9.7) are now 
apphed only for a fixed value of i and a fixed value of j (as opposed to all 1 < z < m and all 

1 < J < ?^), wc obtain that for any integer A the left hand side of (14.9), when regarded as 
a polynomial in g^, admits qi = C,^^ — Xh as a root of multiplicity max(si — (A + 1), 0), and 
as a polynomial in q'j admits q'j = — A/i as a root of multiplicity max(tj — (A + 1), 0). 
Setting h — we obtain the desired divisibility of det M" by the powers of qi — and 
q'. — C^^, and the proof is complete. □ 



15. Random covering surfaces and physical interpretation 

Let qi, ■ ■ ■ ,qn be fixed even integers, and assume 2g = ^^Li Qi ^ 0. Let U he a, fixed 
rhombus on the triangular lattice, centered at the origin and having side I, with I very 
large. Consider the lattice triangular holes of side two W{3ri, 0), . . . , W{3rq, 0), where ri 
is vastly larger than the side of U and is vastly larger than r^-i, for i = 2, . . . ,q. Let 
S be the lattice rhombus rU with opposite sides identified so as to make a torus, where 
r is vastly larger than r^. 

Let C he the collection of all non-overlapping placements C = {Qi, . . . , Qn} of lattice 
triangular holes inside where Qi has side g^, for i = 1, . . . , n. 

Given an integer x > 1, a collection of holes on the triangular lattice T may also be 
viewed as residing on the x-fold refinement of the triangular lattice. Regarded so, 
their correlation is still well-defined; denote it by Ux- 

Let 

^50,= U %{S\QiU---UQnUW{3ri,0)---UW{3rq,0)), 
{Qiv,Q»}eC 

where T^iR) is the set of all lozenge tilings of the region R when regarded as residing on 
the x-fold refinement of the triangular lattice. 
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Then Theorem 1 . 1 imphes that when samphng uniformly at random from Sx , the ratio 
of the probabilities P that Qi, ■ • • , Qn are at mutual distances dij and P' that Qi, • • • , Qn 
are at mutual distances d\j is, in the limit 1 << I « ri « ... « rq « r and for 
large separations between the holes, given by 

2 

p Q-^T,i<i<j<u1i1ii-'^'^dij) 

On the other hand, the Fundamental Theorem of Statistical Physics (see e.g. [9, §40- 
3]) implies that, given a physical system of two dimensional charges of magnitudes 
qiQe, ■ ■ ■ , QnQe (where is the elementary charge), the relative probability that the charges 
are at mutual distances d^j versus d'^^j is given by 

2 

^ ^' (15-2) 

where k is Boltzmann's constant (see e.g. [9, §39-4]), T is absolute temperature, and eo 
is the permitivity of (two dimensional) empty space. 

Thus physical electrostatic interaction in two dimensions is perfectly paralleled by our 
model, with charges corresponding to holes, charge magnitude to the difference between 
the number of right- and left-pointing unit triangles in the hole, and temperature to c/x^, 
where x is the lattice refinement parameter and c is a constant. 

We note that while there is evidence that the superposition principle of Theorem 1.1 is 
independent of the type of the background lattice (see [3, §14]), the specific relationship 
obtained between T and x when doing the above analysis does depend on the lattice. 
We also note that there is an analogous three dimensional model one can consider (see 
[3, §15]), and recent numerical simulations suggests that the parallel to electrostatics also 
holds in three dimensions (see [12]). In the version presented in [3, §15] the relation- 
ship between physical temperature and the lattice refinement parameter is one of direct 
proportionality. 

We conclude this section by providing an equivalent interpretation of our result in 
terms of covering surfaces. 

View the plane with our triangular lattice T as being the (x, j/)-coordinate plane in 
R^. Consider a copy of the lattice L = in placed in such a way that the image of 
the lattice points of L under the orthogonal projection p on the (a;, j/)-coordinate plane 
consists precisely of the lattice points of T. 

Then the triangular lattice lifts under p : L ^ T to the cubic lattice L. Furthermore, 
it is well known that the lozenge tilings of a simply connected region on T lift under p to 
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(15.1) 



Figure 15.1. Lifting with cuts. Figure 15.2. The boundary of a hole becomes a spiral. 

lattice surfaces of L (more precisely, to minimal lattice surfaces — unions S of unit squares 
of L so that the restriction of p to 5" is one-to-one). 

This bijection breaks down if the lattice region has holes. However, it turns out that 
tilings of regions with holes can be lifted to covering surfaces under p. Indeed, consider 
for instance the region with two holes shown in Figure 15.1, and consider a tiling fj, of it 
(also indicated in the figure) . Cut the region along non-intersecting lattice paths leading 
from each hole to the boundary (these cuts are indicated by dotted lines in Figure 15.1). 
The resulting region is simply connected, hence can be lifted to a lattice surface in 
L (four different views of the surface corresponding to the tiling and cuts in Figure 15.1 
are shown in Figure 15.3; the view from the top — with no perspective — looks just like the 
tiling; a further instance is shown in Figure 15.4.). For n G Z, let Sn be the translation of 
Sq by the vector nv, where v is the up-pointing large diagonal of a unit cube of L. Then 
S — UnezSn is independent of the choice of the cuts, and is a covering surface under the 
projection p for the tiling fi. 

The boundary of a hole Q lifts under p : Sn ^ T to a lattice path on L that in general 
is not closed: the vector that separates its endpoints is precisely ch.{Q)v. The full image 
under p~^ of the boundary of the hole — call it a tubular hole of 5" — becomes a collection 
of I ch((5)| identical interspersed spirals of step ch.{Q) (for ch.{Q) — —1 this is illustrated 
in Figure 15.2). This step size is known as the local holonomy of the covering surface. 
The collection of step sizes associated to all the holes form what we call the holonomy 
type of the covering surface. 

By this correspondence, the parallel discussed at the beginning of this section becomes 
a parallel between relative probabilities of two dimensional electrostatic charges being at 
specified mutual distances and relative probabilities of random covering surfaces of given 
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Figure 15.3. Four views of the surface lifting the tihng of Figure 15.1. 



holonomy type (and appropriate periodic boundary conditions) to have their tubular 
holes at specified mutual distances. The parameter x that accounts for temperature 
controls then the "resolution" of the covering surfaces. 

From this perspective, the laws of two dimensional electrostatics arise by averaging 
over all possible discrete geometries of the covering surfaces. 

The joint correlation of the holes is a real number that arises by averaging over all the 
covering surfaces that lift tilings of the complement of the holes. In this paper we studied 
the asymptotic behavior of this number as the holes move away from one another (and 
found that this asymptotic behavior is governed by Coulomb's law for the electrostatic 
potential energy). A complementary question, imperative from the point of view of the 
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Figure 15.4. Two views of the surface lifting a tiling of the 
complement of the holes in Figure 2.1. 



literature on surface models (see e.g. [6], [8] and [7])), is the following. Keep the holes 
fixed, and move around a "probing point," which at each position records the (suitably 
defined) average height of the covering surfaces compatible with the fixed position of the 
holes. What can we say about this average surface in the scaling limit? 

More precisely, let us refine the triangular lattice more and more, keeping the holes of 
fixed side lengths in the refined lattice — this will cause the holes to shrink to points as 
the lattice spacing approaches zero. Then if we arrange for the holes to shrink to points 
chosen beforehand, does the limit of the normalized^ average multivalued height function 
exist? If so, what is this limit? 

^To obtain the normalized height function, as we refine the triangular lattice say ii-fold, we replace 
the lattice in the surface interpretation of tilings presented earlier in this section by the lattice (-^Z)^. 
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Figure 15.5. The limit of the average over ah surfaces of the type shown in Figure 15.4. 

It turns out that the above described hmit surface exists, and is just the horizontal 
coordinate plane with the corresponding points removed. Thus, unlike in the explicit 
limit surface results of Cohn et al. [6] [8] and the related results of [7] [17] [21], to capture 
the interesting behavior we need to consider the un-normalized average height function. 
This occurs also in Kenyon [16], where it is shown that the (un-normalized) average 
height function of Temperleyan domains on the square lattice converges in the scaling 
limit to a harmonic function. 

We show in [5] that the scaling limit of the latter is a sum of "refined helicoids" based 
at the points to which the holes shrink, and with steps proportional to the charge of the 
corresponding holes (if for a helicoid H the points on the same vertical are a distance d 
apart, the s-refinement of H is obtained by taking s — 1 more copies of H, all vertical 
translates of it, and interspersing them so that the points on their union that are on the 
same vertical are a distance d/ s apart). For instance, the scaling limit of the average over 
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all the surfaces of the kind shown in Figure 15.4 is of the type illustrated in Figure 15.5. 

We conclude by mentioning that this limit surface result turns out to be equivalent 
to determining the scaling limit of the discrete vector field defined by the average ori- 
entation of a lozenge covering a fixed unit triangle in the complement of the holes (see 
[5] for details). In this phrasing, the above result states that the scaling limit of this 
discrete field is precisely the two dimensional electrostatic field obtained when regarding 
the holes as charges of magnitude given by the difference between the number of right 
and left monomers in them. Thus the two complementary viewpoints of studying the 
behavior of the correlation when the holes move around, and studying the average height 
function when the holes are fixed, lead to the emergence of two complementary aspects 
of electrostatics: the electrostatic potential energy and the electric field. 
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Appendix. A determinant evaluation 



The object of this section is to prove the following result. 



Theorem A.l. Let si > ■ ■ ■ > Sm > ^ and ti > ■ ■ ■ > tn > 1 be integers. Denote 
S — Si, T — ^^^itj, and assume S > T. Let xi,...,Xm and zi,...,Zn be 

indeterminates. Define N to be the S x S matrix 



N = [A B 



(A.l) 



whose blocks are given by 



0) 



— Xm—^l 

(0 

(-Xm-^l)^ 



(o) 

(-xi-^i)^ 



(si-i) (si-i) 



(-^m-^l)- 



(-a;i-zi)*l+i 



(-^m-^l)*l 

(V) 



('sm+tl-2\ 



(o) 

( — — Zti)^ 



(sl-l) (si-l) 
{-a:i-z„)"l {-a:i-z„)''l+l 



(D 



(?) 



(^") 



(-xi-Zn) 



(-a;i-z„)n+tn-l 



(-Km-^n)*" 

(!l) 



am -1 / 



(A.2) 



and 
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B = 



01 



( ° )x^-'^ ( ' )a 



Then we have 



/S-T-l\ S-T-2 
\ 1 



^S-T-l\„S-T- 



/S-T-l\ S-T-1 
I )-^m 

(S-T-l\ S-T-2 
\ 1 



(A.3) 



detiV 



(A.4) 



Our proof relies on Proposition A. 2 below which presents some interest in its own right. 
Let F be an analytic function, and let f be the ^'-vector 



f=[F{x^),-F'{x,),. 



F{x,),-F'{x2),. 



1 



'(.1-1)! 
1 



F^'^-'\xi), 
F^^^-'\x2), 



]. Tr(sm-i)(j. \] 

[Sm - 1)! 



(A.5) 



where > ■ ■ ■ > 
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For each 1 < i < j < m, consider the following two lists of vectors of type [u, v] , where 
u and V are vectors of lengths Si and sj, respectively. 



The first list, consists of the following Sj groups, totaling Sj{si — Sj + 1) + 1 + 

2 + . . . + (sj — 1) = SiSj — Sj{sj — l)/2 vectors (as before, 0^ denotes the A;- vector with 
all coordinates 0): 



u 



Oo 








• 


■■ 


1 





•• 


• 0^ 


Oo 





G) 


• 


•• 





1 


•• 


• 


Oo 








© • 


•• 








1 • • 


■ ° ' 


Oi 


(S) 





• 


•• 


1 





•• 


• 0^ 


Oi 





© 


• 


•• 





1 


•• 


• 


Oi 










•• 








1 • • 


• > 



Osi—Sj 
Osi—Sj 

Osj — s,- 



( 






1 J 















10 
1 
1 



^ 






> Si 
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0. 
0. 



-8J+1 
-Sj+1 



-Sj + 1 







) 








•• 


1 





■■ 


■ ^ 





•• 





1 


■• 


• 


SJ+3^ . 


•• 








1 • • 


• 



Os,-2 (%"') I 1 ••• 
0s,-2 {'\^) I 1 ■■■ 
o,,_i (%-^) I 1 ••• 0}1. 

(2) 

The second list, L)^j , consists of the foUowing Sj — 1 groups 



(A.6) 



u 




V 
















1 Oo 


-ft) 


1 





•• 


• 


0^ 






1 Oo 


-ft) 





1 


•• 


• 



























1 Oo 


~{sj-l) 








■■ 


• 


1 






1 Oi 


-ft*) 


1 





•• 


• 


0^ 






1 Oi 


-ft*) 





1 


•• 


• 























>Sj-2 




1 Oi 


~(s,-2) 








•• 


• 


1 

J 





Os 

Os. 

0., 



(A.7) 



and has a total of 1 + 2 + . . . + (sj — 1) = Sj{sj — l)/2 vectors. 
Define tfie list to consist of the vectors 

[Os,,...,0^,_,,-u,0^^^,,--- ,Os^_,,v,0^^.^,,--- ,0^^], (A.8) 

where [u, v] ranges first over the list (A. 6) and then over the list (A. 7). Then L^j consists 
of SiSj vectors of length S. 

Let L be the list obtained by concatenating the lists (A.8) as follows: 

L = [Li2, Lis, ■ ■ ■ , L\rni L23, -£'24, ■ ■ ■ , L2 mi • • • 1 L/m—l,m 

]. (A.9) 

Then L consists of ^Kj^^x^ SiSj .S-vectors of type 

[0,...,0,-a,0,...,0,/?,0,...,0]. (A.IO) 

If the two non-zero coordinates of (A.IO) are in positions k and I, k < I, associate with 
the vector (A.IO) the operation 

/, := (A.11) 

Xl - Xk 

on the coordinates of the vector f = [fi, ■ ■ ■■, fs] given by (A. 5). 

Proposition A. 2. After applying the coordinate operations (A. 11) associated to the 
vectors of the list L to the vector f given by (A. 5) and setting Xm = • • • = xi, the 
resulting vector is 

[F(xi), iF'(xi), ^F"{xi), . . . , (A.12) 

Our proof of this proposition is based on the next three Lemmas. They concern the 
special case m = 2, si = a > b = S2 of Proposition A. 2. In this case the vector f becomes 

(A.13) 

Denote by Ra{f, xo){x) the remainder of the truncated Taylor expansion of / around 
Xq: 

fix) = f(xo)+^^{x-xo)+- ■ ■+^^P^^{x-xor-'+{x-xorRa{f,xo){x). (A.14) 
1! a — 1)! 
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Lemma A. 3. Consider the ah coordinate operations (A. 11) corresponding to the vectors 
[— u, v], where the pairs (u, v) are the ones obtained from (A. 6) and (A. 7) by replacing 
Si by a and Sj by b. 

After applying these operations to the vector (A. 13), its a + k + l-st coordinate becomes 



fork = 0,l,...,b-l. 

Proof. We proceed by induction on k. Identify the vectors in the hsts (A. 6) and (A. 7) 
corresponding to Sj = a, Sj = b with the coordinate operations (A. 11) they encode. 

For k = the statement follows from the form of the coordinate operations (A.6)-(A.7) 
and the definition (A. 14) of Ra{F,xi){x2)- 

As indicated in (A. 7), group the b{b — l)/2 operations of (A. 7) in groups of sizes 
6 — 1,6 — 2,..., 1. Note that after the first coordinate operation in the z-th group, 
coordinate a + i+ 1 already reaches its final form. On the other hand, the earliest time 
coordinate a + i + 1 is multiplied by a constant and added to some coordinate that comes 
after it is in the i + 1-st group of (A. 7). Since coordinate a + 1 reaches its final form 
already after operations (A. 6), it follows that by the time coordinate a + i is multiplied 
by a constant and added to row a + k + 1, k > i, row a + i has already reached its final 
form, for alH > 1. 

Therefore, for the induction step we need to prove that 





{X2 - Xi)'' 



1 





{F') + {a)2Ra{F)] 



(■ 



a\ 1 



ij {k-iy. 




+ ■■■ 





Extracting the coefficients of Ra-i(F^'^^) on both sides, (A. 15) is seen to be equivalent 



to 



k-i-l 

E 

3=0 



for i = 0, . . . , /c — 1. After some manipulation (A. 16) turns out to be equivalent to 



n+l 



j=0 



(g + 1 - j)n 
{n+l-j)\j\ 



0, 



(A.17) 



where n = k — i. However, the left hand side of (A.17) can be expressed in terms of 
hypergeometric series as 



n+l 



E(-i)^ 



{a + l-j)n (a + 1). 



{n + l-j)\j\ (n + l)! 



-a, —n — 1 



-a — n 



The resulting 2F1 is zero by Gauss' summation (4.9). 
Lemma A. 4. For /c = 0, 1, . . . , o - 1 one has 



□ 



Ra{F,Xi){x2) 



^-^{i.._.(F(^),.0(^.) 



- {^^{a)^Ra-k+l{F^''-^\x^){x2) + 

+ (-l)'Q(a)fei?a(F,Xi)(x2)|. 



(A.18) 



Proof. By definition (A. 14) we have 



Ra{F, Xi)(x2) = F{X2){X2 - X^)- - J2 ^(^2 - X,)-''+^ . 



i=o 



This implies that 



dx2 



(A.19) 
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Definition (A. 14) also impfies 



Ra-^{F^'\X^){X2) = {^^^^^^2) - F^^\x,) - ^^^^^^{X2 -x^)-■■■ 

1)1^^2 XI) 

for i = 0, . . . , /c. Therefore, writing for brevity x = X2 — Xi, = F«(a;i) and Gi 
F^'^\x2), the right hand side of (A. 18) becomes 



X 



^^ ^k-i Fk_ Fa-i 



+ 



:X 



(a-1)! 



(A.20) 



To prove (A. 18) we need to show that the above expression equals (A. 19). The coeffi- 
cient of Gi in (A.20) is 



k — ij x^~'^ x^ ' 

which clearly agrees with the coefficient of Gi in (A. 19). Agreement of the coefficients of 
Fi on the other hand turns out to amount to 

for i = 0, . . . , /c. 



The left hand side of (A. 21) becomes after some manipulation 

, ,Y^ (-1)^ (a)fc-j ^ (a)fc ^{a)k ^ 



k, i ^ 
-k-{a-l)' 



Gauss' summation (4.9) implies that when a is a large enough negative integer the above 
expression agrees with the right hand side of (A. 21). Then (A. 21) follows since both its 
sides are polynomials in a. □ 
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Lemma A. 5. For all k>0 we have 

—R^{F,xo){x) 



X=Xo 



{a + k)\ 



Proof. Since F is analytic we have by (A. 14) that 

Ra{F, xo){x) = ^ + , . {x - Xo) + 

This imphes 



dx'^ 



Ra{F,XQ){x) 



0!^"^ "^"^ {a + k)\ 1! (a + fc + l)! ^ ' 

Take a; = a;o to obtain the statement of the Lemma. □ 

Proof of Proposition A. 2. We proceed by induction on m. For m = 1 the statement 
is clear. Perform on the vector f the coordinate operations (A. 11) corresponding to 
the initial segment [L12, . . . , Lim] of the list L — a total of S1S2 + ■ ■ ■ + siSm coordinate 
operations. By Lemmas A. 3 and A. 4, the last S — s\ coordinates become 



1 . 1 



1! [Sm - 1)! 

where G{x) = xi)(x). 

By the induction hypothesis, applying the remaining ^2<i<j<m ^i^j coordinate oper- 
ations and letting Xm = ■ • • = X2, the above vector becomes 

[G{X2), ^G'{X2), . . . , Js^G('-''\^^)]- (A-22) 
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However, since G{x) = Rsi{F, xi){x), we obtain by Lemma A. 5 that 



d 



X=Xl 



This shows that when appending the X2 = Xi speciahzation of (A. 22) to the first Si 
coordinates of f, the resulting vector is precisely the one given by (A. 12). □ 

Let F{x) = l/{—x — z) and Gi{x) = x\l> 0. One readily checks that 



1 d!" 
k\ dx^ 



If^'\x) 



l\ 



k + l 
k 



{-X - ^)fc+'+i 



and 



l-k 



for all k,l>0. 

It follows from these equalities that if we let Fj{x) = l/{—x — ^j), j = 0, . . . , n, and 
define the -S-vector h(a;) by 



h{x) = [F,{x), ^Fiix), (^^^^^"'^(^). 



1 



F2{x),-F^ix), 



1 I7(i2-1) 



(t2-i)r 



F„(x),i<(x),. 
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^F^*"-^) (x), Go{x), Gi(x), . . . , Gs-T-i{x)], 

(A.23) 



then the matrix N in the statement of Theorem A.l can 



be written as 



N = 



h{xi) 



h{x2) 



(«2-l) 



L(s^-l)!" {•J^m)^ 

To present our next result, associate to each vector of type (A. 10) in which the non-zero 
coordinates are in positions k < I the row operation 



(A.24) 



Ri 



(3Ri - aRk 

Xi - Xk 



(A.25) 



in matrix N. Then the hst L defined by (A.6)-(A.9) describes Ylii<i<j<m ^i^j oper- 
ations on N. 

In addition, consider the hst 

= [-^127 -^135 • ■ ■ 5 Mn' -^235 -^247 ■ ■ ■ : ■ ■ • ' -^n-l,n]; 

where L'^^ is obtained from the list defined by (A.6)-(A.8) by replacing the m-tuple 
(si, . . . , by the n-tuple (ti, . . . ,tn)- 

Interpret each of the J2i<i<j<n^'i^j vectors of L' as a column operation on matrix N 
using the analog of (A.25) in which the Xj's are replaced by Zi^s. 
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Proposition A. 6. (a). After applying the '^i<i<^j<m ^i^j '^^^ operations (A. 25) of the 
list L to matrix N and setting Xm = ■ ■ ■ = xi = x, the resulting matrix is 



L (S-l)! 



h{x) 
^h'{x) 



(b). After applying the J2i<i<j<n^i^j column operations encoded by L' to matrix N\ 
and setting Zn = • ■ ■ = zi = z, the resulting matrix N2 is of the form 



F{x,z) 

1 dF{x,z) 
1! dx 



1 dF{x,z) 



1! 



dz 



1 d'^F{x,z) 
1!1! dxdz 



1 d'^-^F{x,z) 
(T-1)! 92^-1 

1 d'^F{x,z) 
1!(T-1)! dxdzT-i 



1 * 
1 



d 



L (S-l)! aa^'-i 



'='-^ F{x,z) 1 d'^F{x,z) 

(S-l)!l! dx^-^dz 



1 d^+-^-^Fix,z) p 

(S-1)!(T-1)! dx^-^dz^-i ^ 



where F{x, z) = l/{—x — z). 





(A.26) 



Proof, (a). This follows using the form (A. 24) of matrix N and applying Proposition 
A. 2 to each of its column vectors. 

(b). By (A. 23) the portion of h.^^\x) consisting of its first S — T coordinates has 
the form of the vector f in the statement of Proposition A. 2, for all A; > 0. Therefore, 
applying Proposition A. 2 to each row of A^i we obtain that A^2 has the stated form in the 
portion contained in its first T columns. 

Along the last S — T columns the matrix N2 is the same as A'^i. Since Gi{x) = x^ in 
(A. 23), part (a) implies that the portion of A''! contained in its last S — T columns has 
the form indicated in (A.26). □ 

Lemma A. 7. One has the matrix identity 



i 



0<i,j<n 



I + i 

i + i 
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X 



0<.i,j<n 



3 ^^.i+i+l 



0<i,j<n 



Proof. The stated equality is equivalent to 



n-l 

E 

fe=0 



^ + A [j 



l + kj\k 



l + i + j 



for all < iij < n. This follows by extracting the coefficient of a;* on both sides of the 



identity {x + iy+'{x + ly = {x + iy+'+^ . 

Proof of Theorem A.l. From (8.11) and (8.7) we deduce that 

det N = c = — — 



a=im=i(-^.-S-*^- 



□ 



where c e Q. Write this as 

det TV 



(A.27) 



By Part (b) of Proposition A. 6 we obtain, denoting f{x) = —1/x and using /^"^ (—1) = n\, 
that 



detiV 



ni<i<7<m(^j a^i)*'*^ ni<i<7<n('^J ^i) 



tit. 



x-\ = ■ ■ ■=Xr 



det 



= det 



S-T+k+l 



S-T<k<S-l 
OKKT-1 . 



0<k,l<T-l 



(A.28) 



However, by Lemma A. 7, the last matrix is the product of a lower triangular and an upper 
triangular matrix, both having all diagonal entries equal to 1 — so its determinant is 1. 
Then (A.28) implies that c = 1 in (A.27), thus completing the proof of the Theorem. □ 

Remark A. 8. Theorem A.l generalizes both the Cauchy determinant, which is ob- 
tained when m = n and all Sj's and tj's equal 1, and the Vandermonde determinant, 
obtained when n = 0, Si = • • • = = 1. 

The case m = rn, si = ■ ■ • = Sm = 1, ti = ••• = tn — r is stated as an exercise 
in [18, Ex.42, p.360]. 
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